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I  SUMMARY 


A  computational  shear  band  model  has  been  developed  and  used  to 
predict  the  fragmentation  response  of  explosively-loaded  cylinders  and 
munitions.  To  provide  a  basis  for  the  model,  an  experimental  technique 
was  developed  for  arresting  the  dynamic  expansion  of  exploding  cylinders 
at  desired  loads  (see  Figure  1).  The  technique  was  used  to  infer  the 
phenomenology  and  to  obtain  quantitative  data  for  shear  instabilities. 

Shear  instabilities  were  found  to  nucleate  and  to  grow  in  HF-1  steel  at 
equivalent  plastic  strains  of  about  15%  and  30%,  respectively.  Shear 
instability  surfaces  are  roughly  semicircular,  and  the  surface  displace¬ 
ment  is  proportional  to  the  depth  of  the  shear  band  surface  (see  Figure  14). 

The  computational  model  for  shear  instabilities  was  used  in  conjunction 
with  a  two-dimensional  wave  propagation  code  to  compute  fragment  size 
distributions  from  exploding  cylinders.  The  model  produces  fragments  by 
accounting  for  the  nucleation,  growth,  and  coalescence  of  shear  bands 
according  to  mathematical  expressions  inferred  from  experimental  observa¬ 
tions  and  data.  Shear  surfaces  form  anisotropically  in  the  deforming 
material  according  to  the  three-dimensional  strain  state;  moreover,  as  the 
shear  surfaces  form,  they  relax  the  stresses.  The  resulting  fragment 
size  distributions  (Mott  plots)  are  obtained  by  summing  the  contributions 
from  all  the  computational  cells. 

Four  fragmenting  cylinder  experiments  and  one  fragmenting  projectile 
experiment  differing  in  material,  geometry,  and  explosive,  were  simulated 
computationally  (see  Figures  16  and  17).  Good  agreement  was  obtained 
between  computed  and  observed  fragment  size  distributions  (see  Figures  18, 
20,  and  21).  These  results  indicate  the  basic  correctness  of  the  model 
and  suggest  that  the  model  can  account  for  effects  of  variation  of 
geometry,  strain  rate,  and  pressure.  Further  model  development  should 
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focus  on  early  stages  of  shear  band  nucleation  and  growth,  and  on  the 
stress-strain  relations  of  shear-banded  material. 

Metallographic  examination  of  fragments  from  HF-1  and  tungsten-nickel- 

iron  long-rod  penetrators  showed  that  shear  banding  was  a  principal  mode 

of  projectile  failure  under  oblique  and  normal  impact  conditions  (see 

Figures  22  and  27),  Computational  aimulationa  with  the  two-dimensional 

wave  propagation  code  TROTT  were  made  of  a  plate  impacting  on  edge  a 

o  o  o 

second  plate  at  angles  of  30  ,  45  ,  and  55  .  Both  high  tensile  and  high 
shear  stresses  are  produced  near  the  impacting  edge  of  the  projectile 
plate,  indicating  potential  fracture  and  ahear  banding  regions.  These 
stresses  are  significantly  larger  if  no  slippage  occurs  between  the 
impacting  surfaces.  Future  computations  should  be  performed  with  shear 
banding  and  fracture  subroutines  and  the  results  should  be  compared  with 
metallographic  observations  of  impacted  projectiles. 


II  INTRODUCTION  AND  OBJECTIVES 


Shear  banding,  the  development  of  internal  regions  of  intense 
plastic  strain,  is  a  prominent  mechanism  of  deformation  and  fracture 
leading  to  the  disintegration  of  exploding  rounds,  long-rod  penetrators, 
and  armor  plates.^  The  Ballistics  Research  Laboratory,  the  Navy  Surface 
Weapons  Center,  and  the  Advanced  Research  Projects  Agency  are  therefore 
seeking  to  acquire  an  understanding  of  the  phenomenon,  and  have  contracted 
SRI  to  develop  a  computational  model  that  can  be  used  to  predict  the 
failure  behavior  of  ordnance  structures  and  materials  by  this  mechanism. 

The  model  could  be  used  to  predict  results  of  given  projectile/target 
encounters  and  the  effectiveness  of  fragmenting  rounds  in  given  scenarios, 
thus  allowing  more  accurate  vulnerability  and  lethality  assessments. 
Furthermore,  a  computational  shear  band  model  would  be  useful  in  the  design 
of  improved  weapons  and  armors  and  in  the  development  of  improved  ordnance 
materials. 

SRI  efforts  to  construct  a  shear  band  model  began  in  1974  under 
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contract  with  the  Naval  Surface  Weapons  Center.  The  resulting  model, 
called  SNAG  for  Shear  Nucleation  and  Growth,  treated  the  competing  work¬ 
hardening  and  thermal-softening  effects  of  plastically  deforming  material, 
and  included  plastic  threshold  strain,  shear  band  velocity,  and  critical 
internal  energy  as  parameters.  Reasonable  predictions  of  the  fragmentation 
behavior  of  cylinders  of  HF-1  steel  and  Armco  iron  were  realized.  The  goal 
of  the  research  program  described  here  was  to  refine  this  preliminary  model 
and  apply  it  to  compute  the  fragmentation  behavior  of  exploding  rounds  and 

long-rod  penetrators. 

♦ 

Our  approach  was  to: 

•  Conceive  and  perform  experiments  in  which  shear  bands  were 
produced  and  arrested  at  various  stages  of  their  development. 
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*  Observe  shear  band  phenomenology  and  make  quantitative  measure¬ 
ments  of  such  properties  as  shear  band  numbers,  sizes,  and 
spacings. 

•  Develop  a  computational  model  incorporating  mathematical 
expressions  describing  shear  band  behavior. 

•  Incorporate  the  shear  band  model  in  wave  propagation  codes. 

*  Test  the  model  by  computationally  simulating  experiments  and 
comparing  computed  and  measured  results. 

The  following  section  describes  the  interrupted  exploding  cylinder 
experiments  and  presents  the  resulting  data  on  shear  band  behavior. 
Section  IV  describes  development  of  a  computional  model  based  on  the 
observations  and  data  from  Section  III.  Metallographic  evidence  for 
shear  band  activity  in  long-rod  penetrators  is  presented  in  Section  V, 
and  Section  VI  presents  two-dimensional  slab-slab  impact  calculations 
indicating  the  distribution  of  shear,  bending,  and  tensile  stresses 
expected  in  obliquely  Impacting  rods. 


Ill  SHEAR  BAND  MODEL  DEVELOPMENT  -  EXPERIMENTAL  PHASE 

Our  goal  for  this  phase  of  the  program  was  to  devise  experimental 
techniques  that  would  allow  us  to  study  the  shear  banding  phenomena  from 
nucleation  through  growth  to  coalescence  and  ultimate  fragmentation,  and 
to  quantitatively  measure  parameters  that  would  be  important  in  any 
realistic  predictive  model. 

We  devised  the  "contained  fragmenting  round"  technique  that  enabled 
us  to  load  HF-1  cylinders  at  pressures  and  strain  rates  high  enough  to 
induce  shear  banding,  but  to  contain  the  cylinder  in  such  a  manner  as  to 
arrest  the  shear  banding  process  at  a  level  before  full  fragmentation 
occurs.  We  were  then  able  to  examine  the  material  metallographically , 
observe  the  shear  bands,  and  relate  them  to  the  loading  histories. 

In  the  following  paragraphs,  we  discuss  the  experimental  techniques 
in  greater  detail,  examine  the  results  from  the  experiments,  and  draw 
conclusions  regarding  the  shear  banding  process. 

A.  Contained  Fragmenting  Round  Experimental  Technique 

To  study  the  shear  banding  process  through  its  nucleation-growth- 

coalescence  cycle,  we  needed  a  technique  somewhat  analogous  to  the  plate- 

impact  technique  developed  at  SRI  for  study  of  ductile  and  brittle 
3 

fracture.  The  latter  technique  allows  independent  variation  of  tensile 
stress  magnitude  and  duration,  which  are  the  principal  parameters  con¬ 
trolling  the  nucleation  and  growth  of  ductile  voids  or  brittle  cracks. 

For  shear  banding,  we  expected  that  the  principal  controlling  parameters 
would  be  shear  strain  magnitude  and  duration,  or  alternatively,  shear 
strain  and  strain  rate.  Thus,  we  needed  a  method  that  would  apply  a 
known  shear  load  on  a  specimen  and  then  remove  the  load  after  a  known 
time. 
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We  chose  a  cylindrical  geometry,  as  it  appeared  to  us  to  be  the 
simplest  geometry  for  producing  shear  bands  and  also  the  most  relevant 
to  the  fragmenting  round  problem. 

The  geometry  is  shown  schematically  in  Figure  1.  Dimensions  and 

other  parameters  are  given  in  Table  1.  The  specimen  material  is  HF-1 

steel,  a  material  that  has  been  used  in  naturally  fragmenting  rounds  as 

well  as  in  previous  shear  banding  studies.  The  HF-1  steel  tubes  were 
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provided  by  the  Naval  Surface  Weapons  Center. 

The  specimen  tube  is  filled  with  a  high  explosive  and  detonated 
through  an  H.E.  cone  to  produce  a  relatively  planar  detonation  wave 
through  most  of  the  tube.  The  high  Internal  pressures  from  the  H.E. 
deform  the  tube  outward,  inducing  a  state  of  strain  that  has  a  large 
compressive  radial  component,  a  large  tensile  hoop  component,  and  a 
negligible  axial  component  (ignoring  end  effects).  We  thus  have  a  state 
of  high  shear  strain. 

The  specimen  tube  is  surrounded  by  a  series  of  concentric  containment 
pipes.  The  relatively  thin-walled,  soft  Lucite  buffer  tube  allows  the 
specimen  tube  to  expand  radially  and  prevents  a  sharp  high-pressure  impact 
with  the  main  containment  tube,  which  is  a  very  thick-walled  mild  steel 
pipe.  The  outermost  lead  tube  serves  as  a  momentum  trap;  it  is  precut 
in  the  axial  direction  to  enable  it  to  fly  off  in  several  pieces  in  the 
radial  direction,  taking  most  of  the  radial  momentum  with  it  and  reducing 
the  inward-going  tensile  wave.  Because  of  the  thickness  of  the  containment 
tubes,  the  specimen  tube  experiences  no  radial  tensile  stress  for  at  least 
40  psec,  by  which  time  the  internal  pressures  produced  by  the  H.E.  are 
greatly  reduced,  due  to  escape  of  explosive  gases  from  the  open  ends  of 
the  tube.  The  lack  of  radial  tensile  stresses  inhibits  brittle  fracture, 
leaving  shear  banding  as  the  predominant  failure  mode. 


6 


INITIATION 


CROSS-SECTION  DIAGRAM  OF  CONTAINED  FRAGMENTING  ROUND  EXPERIMENTS 


Table  1 


CONTAINED  FRAGMENTING  ROUND  EXPERIMENTAL  PARAMETERS 


Shot  No. 


5084-2 

5084-3 

5084-4 

5084-5 

Dimensions  (cm  ) 

HF-1  steel  tube 

I.D.  (g) 

7.68 

7.68 

7.68 

7.68 

O.D.  (h) 

12.07 

12.07 

12.07 

1  2.07 

Length  (a) 

IS.  73 

19.05 

19.05 

18.42 

Lucite  sleeve 

thickness  (d) 

1.27 

1.27 

2.86 

1.27 

Steel  containment  pipe 
wall  thickness  (c) 

7.94 

7.94 

6.35 

7.94 

Lead  thickness  (b) 

2.54 

2.54 

2.54 

2.54 

Explosive  column 

Length  (e) 

17.78 

13.97 

14.60 

12.06 

Setback  (f) 

0 

2.54 

2.  22 

2.22 

Explosive  parameters 

H.E,  used 

PETN/  ^ 

Nitro- 

PETN/ 

PETN/ 

p-spheres 

methane 

U-spheres 

u-spheres 

Explosive  density 
(gm/cm3) 

0.98 

1.135 

1.05 

1.34 

Expected  detonation 

Velocity  (mm/psec) 

5.45 

* 

5.75 

6.9 

Pressure  (kbar) 

75 

* 

90 

160 

Note:  (a)-(h)  refer  to 

dimensions  shown 

in  Figure 

1. 

*Data  uncertain  due  to  impurities  in  nltromethane,  but  was  similar  to 
that  in  Shot  No.  5084-2. 

^Pentaerythritol  tetranitrate  (C  H  0  N) /hollow  plastic  micro-spheres 

5  o  JL2  4 

manufactured  by  Emerson  &  Cuming,  Inc.,  Canton,  Mass.;  -  2%  u-spheres 
by  weight. 
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The  strain  and  strain  rate  to  which  the  specimen  tube  is  subjected 
can  be  varied  by  changing  the  explosive  mixture  and  the  thickness  of  the 
Lucite  buffer  tube.  A  higher  pressure  detonation  will  induce  higher 
strain  rates  in  the  specimen  tube,  and  a  thicker  buffer  will  allow  the 
specimen  tube  to  attain  larger  strains.  However,  the  controlling 
parameters  are  not  orthogonal  with  respect  to  the  material  response 
parameters.  For  example,  a  higher  detonation  pressure  'vi.H  produce 
larger  strains  even  if  the  buffer  tube  thickness  remains  the  same. 
Furthermore,  there  is  a  complicated  pattern  of  stress  waves  and  reflec¬ 
tions  while  the  various  containment  tubes  are  "ringing  up"  to  their 
peak  pressures  and  then  unloading.  Hence,  the  loading  history  of  the 
specimen  tube  cannot  be  predicted  easily  without  the  aid  of  computer 
simulations. 

B.  Experimental  Results  -  Shear  Band  Nucleation 

After  running  several  computer  simulations  (using  Cylindrical  PUFF) , 
and  one  proof  test  using  a  4142  steel  specimen  tube  (so  as  not  to  deplete 
our  limited  supply  of  HF-1  steel)  to  establish  the  geometry  of  the 
experiment,  we  performed  four  contained  fragmenting  cylinder  experiments 
with  HF-1  steel,  as  shown  in  Table  1. 

The  first  two  experiments  (Shots  2  and  3)  resulted  primarily  in 
incipient  shear  band  formation,  while  the  latter  two  (Shots  4  and  5) 
resulted  in  substantial  shear  band  growth  and  partial  fragmentation. 

The  specimen  tube  recovered  from  Shot  2  is  shown  in  Figure  2.  The 
partially  fragmented  shear  lip  and  substantial  shear  cracking  at  the 
upper  end  of  the  tube  is  caused  by  the-  large  axial  strain  resulting 
from  the  H.E.  extending  axially  to  the  very  edge  of  the  tube.  Such 
end  effects  were  eliminated  in  subsequent  shots  by  setting  the  H.E. 
back  from  the  tube  edge.  The  tube  recovered  from  shot  3,  shown  in 
Figure  3,  exhibits  a  radial  deformation  that  is  greater  near  the  center 

9 


11 


1 


of  the  tube  than  at  the  enda  because  of  the  explosive  set-back  and 
because  the  internal  pressure  at  the  center  remains  high  for  a  longer 
period  of  tine. 

The  wall  thickness  and  diameter  of  the  recovered  tubes  were  measured 
as  a  function  of  axial  position  and  plotted  to  show  the  final  shapes  as 
shown  in  Figure  4.  The  equivalent  plastic  strain  undergone  by  the  tube 
at  different  axial  positions  was  then  calculated  asstiming  negligible 
axial  strain.  A  definition  of  equivalent  plastic  strain  and  details  of 
this  calculation  can  be  found  in  Appendix  A.  Equivalent  plastic  strain 
was  selected  as  a  parameter  relevant  to  the  shear  banding  process  that 
could  be  compared  with  computer  simulations.  For  Shot  3.  the  strain  at 
the  position  of  maximum  deformation  was  «  28%. 

Other  than  the  edge  effect  damage  discussed  above,  the  only  visible 
sign  of  material  failure  in  Shots  2  and  3  was  the  presence  of  minute 
axial  striations  on  the  inner  tube  surface,  as  can  be  seen  in  Figures  2(b) 
and  3(b).  These  striations  appear  after  the  tubes  have  deformed  to  an 
equivalent  plastic  strain  of  15%  and  often  extend  for  an  inch  or  more 
along  the  cylinder  axis.  We  hypothesized  that  these  striations  are  the 
sites  for  subsequent  shear  banding.  To  test  the  hypothesis,  we  made  a 
cut  in  the  tube  recovered  from  Shot  3,  perpendicular  to  the  cylinder 
axis  near  the  point  of  maximum  deformation.  Photomicrographs  of  the 
polished  and  etched  surface  are  shown  in  Figure  5,  and  these  pictures 
clearly  show  the  incipient  stages  of  shear  band  formation.  The  first 
stage  appears  as  surface  striations,  as  the  material  experiences  non- 
uniform  radial  deformation.  Further  radial  expansion  causes  shear 
motion  to  concentrate  at  the  striations.  Once  this  nonuniform  shear 
deformation  has  begun,  the  local  material  becomes  hotter  and  weaker, 

o 

and  the  shear  band  continues  to  propagate  into  the  material  along  the  45 
plane  of  maximum  shear  stress,  provided  that  the  internal  pressures 
remain  high  enough  to  continue  radial  expansion.  The  material  along  the 
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FIGURE  4  OUTLINES  OF  HF-1  STEEL  TUBE  WALLS  BEFORE  AND  AFTER  DEFORMATION 
IN  CONTAINED  FRAGMENTING  ROUND  EXPERIMENTS 


0.5  mm 


MP -5084- 15 


FIGURE  5  INCIPIENT  SHEAR  BAND  FORMATION  IN  HF-1  STEEL  SHOWING 

(a)  SURFACE  DIMPLING,  (b)  INITIAL  SHEAR  MOTION  AT  A  DIMPLING 
SITE.  AND  (cl  SUBSEQUENT  GROWTH 
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shear  band  has  a  reduced  tensile  strength,  so  that  subsequent  tensile 
stresses  can  produce  a  crack  in  the  shear-banded  zone,  as  seen  in  Figure 
5(c). 

Other  cuts  were  made  in  the  HF-1  steel  tubes  recovered  from  Shot  3 
and  from, the  regions  of  Shot  2  unaffected  by  edge  effect  damage,  but 
metallographic  observation  revealed  no  damage  other  than  the  type  shown 
in  Figure  5--incipient  shear  band  formation  and  band  growth  to  a  depth 
of  only  a  few  millimeters.  Thus  we  can  conclude  that  in  this  geometry 
and  at  the  strain  rates  produced  by  H.E.  detonation  pressures  on  the 
order  of  100  kbar,  shear  bands  begin  to  nucleate  at  equivalent  plastic 
strains  greater  than  ss  15%,  but  have  not  grown  to  significant  size  (the 
size  where  coalescence  or  fragmentation  may  take  place)  at  equivalent 
plastic  strains  of  28%. 
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C.  Experimental  Rwulti  -  Shear  Band  Orientation,  Growth,  and  Coalescence 

Contained  fragmenting  round  experiments  4  and  5  resulted  in  much 

greater  levels  of  damage  to  the  specimen  tube  than  the  earlier  shots. 

This  was  caused  by  higher  detonation  pressures  of  the  H.E.  (Shot  5)  or 

by  thicker  Lucite  buffer  tubes  (Shot  4).  In  both  cases  the  containment 

system  was  not  strong  enough  to  keep  the  specimen  tube  intact.  As  a 

result,  the  HF-1  steel  fragmented  into  many  relatively  large  pieces, 

most  of  which  were  recovered.  Figure  6  is  a  photograph  of  recovered 

fragments  from  Shots  4  and  5;  recovered  fragment  size  distributions  are 

shown  in  Figure  7.  (The  equivalent  radius  (R  )  is  the  radius  of  a  sphere 

£ 

having  the  same  mass  as  the  fragment.) 

Figure  8  depicts  cylindrical  tube  elements  containing  shear  bands 
of  the  three  possible  orientation  types  for  cylindrical  symmetry.  For 
every  shear  band,  the  material  on  one  side  of  the  band  slips  with  respect 
to  the  material  on  the  other  side.  If  a  shear  band  intersects  a  material 
surface,  part  of  the  slip  plane  is  exposed  above  the  surface  (the  shaded 
regions  in  Figure  8).  We  define  t  as  the  length  of  this  exposed  region, 
or,  more  generally,  as  the  length  of  the  band  in  the  direction  perpen¬ 
dicular  to  the  slip  direction.  If  we  cut  the  material  along  a  plane 
perpendicular  to  the  shear  band  length,  we  see  the  band  edge-on  and  can 
define  the  parameters  d  as  the  depth  of  the  band  in  the  slip  direction, 
or  the  average  depth  of  the  material  on  opposite  sides  of  the  band,  and 
B  as  the  shear  displacement  or  slip  magnitude  at  the  surface.  (B  cor¬ 
responds  to  Burger's  vector  in  atomic  dislocation  theory.) 

For  orientation  type  1,  l  is  in  the  axial  direction; 

for  type  2,  l  is  in  the  circumferential  direction;  and  for  type  3, 

l  la  in  the  radial  direction.  Within  each  orientation  type,  there  are 

o 

two  possible  directions — 90  with  respect  to  each  other — representing 
planes  of  shear  strain  symmetry.  All  six  of  these  directions  represent 
slip  along  planes  of  maximum  shear  strain.  For  a  radially  expanding 
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FIGURE  6  HF-1  STEEL  FRAGMENTS  RECOVERED  FROM  (a)  SHOT  5084-4 
AND  (b)  SHOT  5084-5 

Arrow  shows  fragment  used  for  quantitative  shear  band  study. 
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FIGURE  7  FRAGMENT  SIZE  DISTRIBUTIONS  FROM  (a)  SHOT  5084-4  AND  (b)  SHOT  5084-5 


ORIENTATION  TYPE  1 


8  -  Length  in  Direction 
Perpendicular  to 
Slip  Motion 

d  —  Depth  Along  45° 
Slip  Plane 

B  -  Shear  Displacement 
Along  Slip  Plane 

Z  —  Axial  Direction 

r  —  Radial  Direction 

8  —  Circumferential 
Direction 


MA-5084-21 

FIGURE  8  GEOMETRY  AND  NOMENCLATURE  FOR  SHEAR  BANDS  IN  CYLINDRICAL 
ROUND  EXPERIMENTS 


cylindrical  tube,  the  largeet  of  theae  ahear  strains  are  in  the  plane 
perpendicular  to  the  axial  direction,  so  shear  bands  of  orientation  type 
1  are  predominant  in  the  fragmenting  cylinder  experiments.  Type  1 
accounts  for  all  of  the  shear  band  damage  in  Shot  4,  as  well  as  most 
of  the  damage  in  Shot  5. 

* 

Figure  9  shows  a  band  that,  for  much  of  its  length,  has  not 
separated  in  tension,  although  the  material  on  opposite  sides  of  the 
transformed  band  has  slipped  distances  that  are  orders  of  magnitude 
greater  than  the  5-p.  band  thickness.  Large  displacements,  therefore, 
do  not  in  themselves  cause  the  shear  band  to  crack.  Figure  10,  on  the 
other  hand,  shows  a  band  that  has  almost  completely  separated.  This 
separation  occurs  sometimes  entirely  within  the  transformed  band,  some¬ 
time  adjacent  to  the  band  but  parallel  to  it,  and  sometimes  across  the 
band  at  a  slight  angle.  Whether  a  crack  lies  entirely  within  a  shear 
band  is  probably  determined  by  the  time  elapsed  between  band  formation 
and  subsequent  tensile  failure.  If  the  hot  shear  band  has  had  time  to 
be  quenched  by  the  surrounding  material  and  thus  regain  some  of  its 
strength,  the  crack  will  more  likely  lie  partially  outside  the  band.  In 
any  case,  shear  cracking  is  seen  only  in  the  vicinity  of  shear  banding 
in  the  fragmenting  cylinder  experiments,  and  so  for  the  quantitative 
damage  assessment  discussed  below,  shear  bands  and  shear  cracks  are 
considered  the  same. 

The  fragments  recovered  from  Shot  4  were  created  almost  entirely  by 
the  shear  banding  process;  the  lateral  edges  of  the  fragments  were  all 
ahear  bands  of  orientation  type  1.  The  fragments  recovered  from  Shot  5, 
on  the  other  hand,  showed  shear  band  damage  only  in  the  region  from  the 
inner  surface  to  about  halfway  through  the  wall  thickness,  and  showed 

* 

The  band  appears  white  after  etching  because  the  steel  has  been  metal- 
lographically  transformed  by  the  heat  generated  in  the  band  formation 
process. 
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INSIDE  CYLINDER  SURFACE 


FIGURE  9  DETAILED  PHOTOMICROGRAPH  OF  CRACK  9A  ON  CUT  7,  WITH  INSET  SHOWING 
HIGH  MAGNIFICATION  OF  SHEAR  BAND  REGION 


INSIDE  CYLINDER  SURFACE 


FIGURE  10  DETAILED  PHOTOMICROGRAPH  OF  CRACK  3A  ON  CUT  4,  SHOWING  VARIOUS 
RELATIONSHIPS  BETWEEN  SHEAR  BANDING  AND  SUBSEQUENT  CRACKING 


primarily  brittle  fracture  the  remainder  of  the  way  toward  the  outer 
* 

surface.  We  therefore  chose  to  study  in  detail  a  fragment  from  Shot  4 
and  selected  the  one  denoted  by  an  arrow  in  Figure  6,  as  it  appeared  to 
contain  a  large  number  of  shear  bands  that  intersected  the  inner  surface. 

Photographs  of  the  fragment  from  two  different  directions  are  seen 

in  Figure  11.  The  two  directions  A  and  B  highlight  the  shear  bands  that 

o 

intersect  the  inner  surface  at  the  two  different  45  planes.  For  identi¬ 
fication  purposes,  the  shear  bands  are  labelled  by  a  number  and  a  letter 
(e.g. ,  band  IB  is  the  band  numbered  1  that  is  visible  from  direction  B) . 
To  allow  examination  of  internal  damage,  the  fragment  was  cut  into  thin 
layers  perpendicular  to  the  axial  direction  along  the  lines  shown  in 
Figure  11.  The  cut  surfaces  were  then  polished,  etched,  and  examined 
metallographically.  Figure  12  is  a  composite  photomicrograph  of  the 
surface  of  cut  4,  showing  clearly  how  the  shear  bands  visible  on  the 
inner  surface  extend  into  the  material,  and  also  revealing  internal  shear 
bands. 

Examination  of  the  composite  photomicrographs  taken  from  each  of  the 
eight  cuts  revealed  that  more  than  two-thirds  of  the  shear  bands 
intersected  the  inner  surface  of  the  cylinder.  Two  bands  extended  from 
the  interior  to  the  outer  surface,  and  the  remainder  lay  entirely  within 
the  interior.  We  believe  that  the  bands  originated  at  their  innermost 
edge  where  the  shear  strains  were  larger  and  occurred  earlier,  and  pro¬ 
ceeded  to  grow  outward;  however,  we  cannot  be  certain  of  that  in  each 
case.  Although  most  of  the  bands  in  this  fragment  stopped  without  inter¬ 
acting  with  any  other  bands,  some  bands  did  intersect  others.  Some  short 
bands  intersected  and  stopped  at  the  middle  of  a  longer  band,  indicating 
that  the  short  band  was  formed  at  a  later  time  than  the  longer  band. 

* 

This  was  probably  caused  by  the  premature  failure  of  the  containment 
system,  which  allowed  radial  relief  waves  to  propagate  into  the  HF-1 
tube,  causing  tensile  failures. 
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FIGURE  11  HF-1  STEEL  FRAGMENT  FROM  SHOT  5084-4  SELECTED  FOR  USE  IN 

QUANTITATIVE  SHEAR  BAND  STUDY,  PHOTOGRAPHED  FROM  (a)  DIRECTION 
A  AND  (b)  DIRECTION  B 


FIGURE  12  COMPOSITE  PHOTOMICROGRAPHS  OF  CUT  4,  SHOWING  VARIOUS  POSITIONS  AND 
ORIENTATIONS  OF  SHEAR  BANDS  AND  CRACKS  WITHIN  HF-1  FRAGMENT 

Note  shear  bands  14A  and  7B  intersect  to  form  tensile  crack. 


Other  intersecting  bands,  such  as  bands  14A  and  7B  in  Figure  12  appeared 
to  grow  simultaneously  and  interact  to  form  a  tensile  crack  beginning 
at  their  point  of  intersection. 

The  depth  (d)  of  all  of  the  bands  revealed  on  the  eight  composite 

photomicrographs  was  measured,  yielding  the  distribution  shown  in 

Figure  13.  Since  the  cuts  did  not  necessarily  intersect  the  bands  at 

their  point  of  maximum  depth,  the  depth  values  measured  are  only  lower 

bounds.  If  we  examined  a  greater  number  of  bands,  a  statistical  trans- 

4 

formation  similar  to  that  used  for  tensile  failure  analysis  could  be 
used  to  determine  the  actual  distribution. 

The  surface  shear  displacement  (B)  for  all  the  bands  intersecting 
the  surface  was  measured  and  the  relationship  between  d  and  B  for  all  of 
those  bands  is  plotted  in  Figure  14.  Three  bands  that  extended  over 
several  cuts  are  identified  by  special  symbols.  Most  of  the  points  lie 
close  to  the  line  d  =  15B,  indicating  a  direct  proportion  between  the 
shear  band  depth  and  surface  displacement. 

The  shear  band  depths  were  compared  with  the  length  at  the  inner 
surface  so  as  to  determine  the  approximate  shape  of  the  shear  band  in 
its  slip  plane.  The  results  are  plotted  in  Figure  15,  along  with  the 
locus  of  points  for  a  semicircular  shear  band  shape.  Again,  noting  that 
the  depth  values  are  only  lower  bounds,  it  appears  that  l  is  on  the 
order  of  twice  d.  If  a  shear  band  is  assumed  to  start  growing  at  a 
single  point,  then  it  appears  to  grow  in  depth  about  as  much  as  it  grows 
in  surface  length  on  each  side  of  its  starting  point.  Without  an 
extensive  program  of  successive  cycles  of  grinding,  polishing,  etching, 
and  photography,  we  cannot  observe  the  precise  shape  of  the  outer  edges 
of  the  shear  bands.  However,  since  d  *=»  15B,  we  can  get  an  approximate 
idea  by  observing  the  shape  of  the  slip  plane  exposed  above  the  inner 
surface  in  Figure  11  (or  in  other  words,  how  B  varies  along  the  length). 
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NUMBER/cm3  SHOWING  DEPTH  GREATER  THAN  THE  ABCISSA 


MA-5084-1 2 

FIGURE  13  DEPTH  DISTRIBUTIONS  OF  SHEAR  BANDS  IN  HF-1  STEEL  FRAGMENT 
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SHEAR  BAND  DEPTH  ALONG  45°  SLIP  PLANE  (cm) 


MAXIMUM  DISPLACEMENT  OF  SHEAR  BAND 
ALONG  SLIP  PLANE  (cm) 

MA-5084-10 


FIGURE  14  SHEAR  BAND  DEPTH  (d)  VERSUS  DISPLACEMENT  (B)  DATA 
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LENGTH  PARALLEL  TO  CYLINDER  AXIS  (cm) 


SHEAR  BAND  DEPTH  ALONG  SLIP  PLANE  (cm) 

MA-5084-1 1 

FIGURE  15  SHEAR  BAND  LENGTH  (8)  VERSUS  DEPTH  (d)  DATA 
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Using  this  method,  we  find  that  some  shear  bands  appear  semicircular, 
others  seem  elliptical,  and  others  have  more  complex  shapes.  But  for 
computational  modeling,  it  appears  that  the  bands  can  be  represented  as 
semicircles. 

D.  Conclusions 

The  above  work  has  lead  to  the  following  experimental  observations: 

•  Multiple  waves  of  tension  and  compression  propagate  through 
the  thickness  of  a  fragmenting  cylinder. 

•  Shear  bands  nucleate  and  grow  in  proportion  to  the  amount  of 
shear  strain. 

•  Bands  usually  appear  to  start  at  the  inner  diameter  of  fragmenting 
rounds,  but  often  they  appear  at  some  depth  within  the  material. 

•  For  bands  with  large  displacement,  the  radial  and  axial  extent 
of  the  bands  are  comparable;  hence,  bands  are  approximately 
circular. 

•  The  maximum  dislocation  of  material  across  a  band  occurs  at 
the  center  and  is  5  to  10%  of  the  radial  or  axial  extent  of  the 
band. 

A  o 

■  Bands  form  at  45  to  the  principal  directions;  that  is,  in  the 
directions  of  maximum  shear  stress, 

•  During  plastic  flow,  many  small  ripples  or  striations  appear 

on  the  inner  surface  at  a  strain  of  about  15%.  At  larger  strains, 
a  small  number  of  these  striations  form  shear  bands. 

•  Brittle  cracks  often  form  along  the  bands  at  a  late  stage. 

In  conclusion,  the  contained  fragmenting  round  technique  has  proven 
successful  in  increasing  our  knowledge  of  shear  banding  phenomena.  How¬ 
ever  more  experiments  are  required  to  better  establish  the  nucleation 
criteria  and  the  band  growth  velocities. 
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IV  SHEAR  BAND  MODEL  DEVELOPMENT  -  COMPUTATIONAL  PHASE 


The  computational  model  for  shear  banding  is  a  detailed  description 
of  the  processes  by  which  shear  bands  form  and  grow  in  a  material  under¬ 
going  large  deformation.  This  section  contains  a  review  of  the  phenomena, 
an  analysis  relating  the  nucleation  and  growth  to  material  properties,  and 
a  discussion  of  additional  features  needed  to  implement  the  model  in  a 
subroutine.  A  listing  of  the  subroutine  and  further  information  about 
the  parameters  are  given  in  Appendix  B. 

The  subroutine  containing  shear  banding  acts  as  a  stress-strain 
relation  in  a  one-  or  two-dimensional  wave  propagation  computer  program. 
The  program  computes  strain  increments  in  a  cell  of  material  and  then 
calls  the  shear  band  subroutine  to  calculate  the  corresponding  stresses. 

If  a  threshold  criterion  has  been  exceeded,  the  subroutine  also  calculates 
the  growth  of  existing  bands  and  the  nucleation  of  new  bands,  and  checks 
whether  coalescence  and  fragmentation  are  occurring.  For  each  cell  of 
material,  there  are  data  arrays  containing  the  size,  number,  and 
orientation  of  all  bands. 

The  shear  band  model  is  intended  to  account  for  the  experimental 
observations  listed  at  the  end  of  the  previous  section.  In  addition  to 
those  features,  the  model  includes  rate-dependent  elastic-plastic  stress- 
strain  relations. 

A.  Derivation  of  the  Model 

The  model  includes  nucleation  of  shear  bands,  growth  of  the  bands, 
coalescence  to  form  fragments,  and  relaxation  of  the  stresses  in  accord¬ 
ance  with  the  growth  of  damage.  In  the  following  discussion,  the  physical 
processes  of  failure  by  shear  banding  are  outlined.  Then  each  process 
is  examined  and  equations  are  derived  to  represent  the  process  in  the 
model. 
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Under  a  high  shear  strain,  materials  first  deform  uniformly, 
following  approximately  the  usual  assumptions  of  homogeneity  and  isotropy. 
But  as  the  plastic  strains  become  large,  inhomogeneities  may  develop.  At 
these  inhomogeneities,  a  slip  like  a  small  crack  can  occur  parallel  to 
the  direction  of  maximum  shear  stress  and  strain.  With  continued 
straining,  the  slipped  region  or  shear  band  grows  in  a  circular  fashion. 
That  is,  the  growth  rate  in  the  direction  of  the  shear  slip  is  approxi¬ 
mately  equal  to  that  in  the  direction  perpendicular  to  the  slip.  Following 
nucleation  of  the  bands,  growth  can  occur  until  the  bands  coalesce.  On 
a  cross  section  the  band  is  observed  as  a  narrow  strip  of  highly  deformed 
material.  Probably  the  strip  is  heated  approximately  to  melting  during 
growth  of  the  band.  At  the  end  of  growth,  the  band  material  remains 
welded  together,  but  with  a  reduced  tensile  strength.  Hence,  subsequent 
tensile  stresses  in  the  material  will  fracture  in  or  along  the  band, 

1.  Nucleation 

The  nucleation  process  includes  several  features:  a  threshold 
for  nucleation,  the  rate  of  nucleation,  the  size  of  bands  at  nucleation, 
and  the  orientation  of  the  bands.  When  the  plastic  shear  strain  exceeds 
a  threshold  value,  nucleation  begins.  In  the  HF-1  material  there  is  a 
threshold  at  about  15%  at  which  small  irregularities  appeared  in  the 
surface.  But  since  significant  slipping  begins  at  30%,  this  value  of 
plastic  strain  serves  as  the  threshold  criterion.  Thus,  the  threshold 
is  an  observed  quantity  and  probably  varies  from  one  material  to  the 
next. 

2 

The  derivation  conducted  by  Seaman  et  al.  indicated  that  the 
rate  should  be  proportional  to  the  square  of  the  yield  strength  and  to 
the  applied  strain  rate;  the  resulting  calculations  correlated  fairly 
well  with  the  measured  fragment  size  distributions.  In  that  study  the 
orientation  of  the  shear  bands  was  not  considered,  so  nucleation  referred 
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to  the  sun  of  the  bands  In  all  directions  and  the  shear  strain  was  a  scalar 
(undirected)  quantity. 

Here  a  nucleation  rate  function  is  proposed  to  account  for  the 
observed  relationships  to  plastic  shear  strain  and  yield  strength  and 
to  provide  an  orientation  dependence  for  the  shear  band  process.  The 
major  assumption  in  the  following  proposed  expression  is  that  the  nuclea¬ 
tion  rate  in  a  direction  is  proportional  to  the  plastic  shear  strain  in 
the  same  direction. 


dt 


C  F 
n  n 


PS 

de  cp6 
dt 


(1) 


p 

Here  E  is  a  critical  energy,  taken  to  be  the  melt  energy;  dE  dt 
cr 

is  the  rate  of  increase  of  plastic  strain  energy;  is  the  number  of 

.  *  ps 

shear  bands  per  unit  volume  in  the  Cp8  direction;  and  is  the 

plastic  shear  strain  in  the  same  direction.  C  is  a  coefficient  with 

n 

the  units  of  time  squared,  and  F  is  a  factor  representing  the  fraction 

n 

of  solid  angle  assigned  to  the  orientation.  Orientation  effects  are 

p 

treated  in  more  detail  later.  The  factor  dE  /dt  is  proportional  to  the 
product  of  the  shear  strain  rate  and  the  yield  strength;  hence,  the 
nucleation  rate  is  proportional  to  the  cube  of  the  strain  rate. 

The  plastic  strain  energy  is  calculated  as  follows. 

dEp  =  V  E  <x'  dfP  (2) 

X  J  i  J 

where  V  is  specific  volume, 

*  . 

Here  cp  i*  the  polar  angle  and  8  is  the  angle  of  rotation  about  the 
pole. 
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is  the  deviator  stress  in  the  ij  direction,  and 

p 

is  the  plastic  strain  increment  in  the  ij  direction. 


The  plastic  strains  in  any  orientation  are  given  by  the  flow  rule 


”rith  the  ^/l.5  factor  in  the  definition  of  d€P,  dfP  is  the  shear  strain 
in  any  geometry  where  two  strain  components  are  equal.  The  plastic 
shear  strains  in  the  six  orientations,  dfP®  ,  are  obtained  by  using 
f°r  Ojj  in  Eq.  3  only  the  shearing  component  of  the  deviator  stress 
acting  on  each  of  the  six  planes.  These  shear  stress  components  are 


listed  in  Table  2 


Table  2 


SHEAR  STRESS  COMPONENTS  ACTING  ON  SHEAR  BAND  PLANES 


Direction 

Index 


Direction 


S.B.  Type 


Shear  Stress  Component 


1 

axial 

— 

a' 

12 

2 

radial 

- 

/ 

a 

12 

3 

circumferential 

- 

none 

4 

o 

45  between 

axial  and  radial 

2 

-  a'  ) 

2  11  22 

5 

o 

45  between  axial 

and  circumferential 

3 

-  °33>2 

+  2 

d'2 

12 

6 

o 

45  between  radial 

and  circumferential 

1 

l  -  °33>2 

+  2 

a'2 

12 

Notes 

a  The  direction  index  indicates  one  of  6  stress  directions, 
b  Shear  band  orientation  type  as  shown  in  Figure  8. 

c  Orientations  used  in  a  two-dimensional,  axisymmetric  code  (TROTT) . 
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Shear  bands  may  appear  in  many  orientations  and  with  many  sizes. 


However,  for  one-  and  two-dimensional  problems,  we  have  chosen  six 

specific  orientations  to  represent  the  continuum  of  all  orientations. 

The  orientations  provided  are  for  bands  normal  to  three  perpendicular 

o 

coordinate  directions  and  for  three  more  at  45  to  the  first  three.  The 
latter  three  lie  in  the  x-y  plane,  x-z  plane,  and  y-z  plane,  respectively. 
For  each  of  these  orientations  a  fraction  of  all  possible  orientations 
is  assigned.  Each  of  the  first  three  orientations  represents  two  points 
on  the  unit  sphere  (for  example,  x  =  ±  1)  whereas  each  of  the  45°  orienta¬ 
tions  represents  four  points.  Therefore,  it  seemed  reasonable  to  weigh 

the  six  points  with  factors,  F  =  1/9,  1/9,  1/9,  2/9,  2/9,  and  2/9, 

n 

respectively. 

For  each  orientation  there  is  a  set  of  size  groups  consisting 
of  a  number  of  bands  per  cubic  centimeter  with  a  radius  .  Initial 

ly  these  size  groups  represent  the  continuous  distribution 


N  =  N  exp(-R/R  )  (6) 

go  n 

3 

where  N  is  the  number  of  bands/cm  greater  than  R 

g 

N  is  the  total  number 
o 

R  is  a  nucleation  size  parameter, 
n 

This  continuous  distribution  is  divided  into  size  groups  such  that  the 
th  . 

i  group  contains  AN^  bands  given  by  the  following  formula: 


An  =  fTAt  [exp(-R  /R  )  -  exp(-Ra  ,/R  )] 
i  in  i+1  n 


(7) 


By  using  size  groups  that  are  discrete  in  size  and  orientation  it  is 
possible  to  represent  growth  and  nucleation  processes  that  are  strongly 
dependent  on  shear  band  size  and  plastic  strain  tensors  that  vary  greatly 


with  orientation 
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2.  Growth 

The  growth  of  a  shear  band  refers  to  the  increase  in  radius  of 

2 

the  slipped  area.  Following  Seaman,  et  al.,  we  have  assumed  that  the 
velocity  of  the  shear  band  is  proportional  both  to  its  radius  and  to  the 
plastic  shear  strain  rate;  that  is, 


dRi  i 

dt  “  CG 


where  C  is  a  growth  coefficient,  d€ ^o/dt  is  the  plastic  shear  strain 
G 

t  h 

rate  on  the  cpO  plane,  and  R.  is  the  radius  of  the  i  band  in  a 
size  distribution.  The  growth  in  Eq.  (8)  is  limited  by  two  additional 
requirements:  The  velocity  cannot  exceed  the  shear  wave  velocity,  and 
the  total  distortion  of  all  the  bands  in  the  cp0  orientation  cannot 
exceed  the  total  plastic  shear  strain  on  the  orientation.  The  first  of 
these  requirements  is  simply  stated  as 


R„  s  R  +  V  At 
2  l  max 


where  V  is  the  shear  wave  velocity.  The  second  requirement  is  based 
max 

on  a  restatement  of  the  Orowan  equation  used  in  dislocation  dynamics. 

The  appropriate  expression  here  is 


—  (ff  S  N  R2B  ) 

dt  r  i  i  i 


where  is  the  average  shear  slip  over  the  area  ff  R^  of  the  band. 

If  (the  number  of  shear  bands)  is  replaced  by  the  band  edge  length 
=  2ffR^N^,  and  is  presumed  constant  during  differentiation,  the 

usual  Orowan  equation  is  obtained.  Over  a  finite  time  interval  and  for 


a  fixed  number  of  bands,  Eq.  9  takes  the  form 


*P  -  «  VVL  -  B11RL> 


where  subscripts  1  and  2  indicate  quantities  before  and  after  the  time 
step.  If  AC  from  Eq.  10  exceeds  the  total  imposed  increment  in 
plastic  strain,  then  R  is  reduced  to  meet  the  requirement. 

A 

From  the  experimental  data  of  our  fragmentation  tests  4  and 
5,  the  dislocation  in  a  circular  shear  band  appears  to  be  5  to  10%  of 
the  radius.  That  is, 


B  =  bR  . 


Thus,  the  dislocation  increases  in  proportion  to  the  radius.  Eq.  10 
then  takes  the  form 


3.  Effects  of  Damage 

Damage  caused  by  shear  bands  leads  to  a  relaxation  of  the 
stresses  as  well  as  to  eventual  separation  of  the  material  into  fragments. 
The  relaxation  effect  is  produced  when  damage  reduces  the  area  across 
which  stress  can  be  transmitted.  The  average  stress  cr  over  the  total  or 
gross  section  is  then  related  to  the  stress  (Ts  in  the  intact  or  solid 
material  through  a  damage  quantity  D.  In  tension  the  relation  is 


a  =  (Xs  (1- D)  . 


The  gross  section  is  composed  of  a  solid  area  and  a  damaged  area 
Ag(l-D) .  The  next  step  is  to  derive  an  appropriate  expression  for  D 
that  will  approach  1.0  for  full  fragmentation  and  account  for  the 
anisotropy  of  damage. 


i 


To  determine  an  appropriate  damage  function,  we  examine  two 
quantities:  the  total  plastic  strain  taken  by  the  bands  and  the  frac¬ 

tion  of  the  material  that  has  fragmented.  The  total  plastic  strain  in 
any  orientation  is,  according  to  Eq.  12, 

€P  =  ffb  p  (14) 

When  the  material  breaks  into  fragments,  the  volume  of  fragments  is 

V  =  T  £  N  f(nf>3  =  1.0  (15) 

F  i  i 

where  T  is  a  dimensionless  volume  factor  (about  4.0) 

F 

f  f  th 

N^  and  R^  are  number  and  radius  of  fragments  in  the  i  size 

group . 

The  numbers  of  fragments  are  related  to  the  numbers  of  bands  through  the 
factor  P: 

N*  =  P  N  .  (16) 

Chunky  fragments  usually  have  six  or  eight  sides  that  were  each  formed 
by  a  band  or  crack.  Because  each  crack  forms  a  side  of  two  fragments, 
three  or  four  cracks  must  be  associated  with  each  fragment.  Therefore, 

P  equals  1/3  or  1/4.  Similarly,  the  fragment  sizes  are  related  to  the 
band  sizes  through  the  factor  y-. 


R 


f 

i 


=  y  R 


(17) 


f  f  3 

Here  R  is  defined  so  that  the  fragment  volume  is  T  (R  )  where  T 
i  F  i  F 

is  about  41T/3.  The  cracks  forming  the  fragment  sides  have  about  the 
same  area  as  the  sides,  so  y  is  approximately  equal  to  1.  Now  Eq.  15 
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can  be  rewritten  as  a  double  sum  of  shear  bands  over  orientations  and 
size  groups: 


A  comparison  of  Eqs.  14  and  18  shows  that  damage  appears  to  be  related 
3 

to  NR  at  all  times  during  the  calculation.  We  will  use  Eq.  18  as  our 

definition  of  complete  fragmentation  (V  =  1)  and  let  V  be  the  fragmented 

fraction  for  values  less  than  full  fragmentation.  Then  the  quantity 
3  3 

TfP7  £  NCp01RCpQi  is  a  measure  of  the  damage  associated  with  the  bands 

in  the  orientation.  We  can  assume  that  the  damage  felt  in  any 

2 

direction  is  proportional  to  cos  Ip  where  Ip  is  the  angle  between  the 
stress  direction  and  the  damage  orientation.  Then  the  damage  in  the 
x  direction  is 


D 

xx 


=  3T  Jy[T 
F  xx 


+ 


0.5(T 


xy 


+ 


T  )] 


xz 


(19) 


where 


=  3iU  E  ' 
i  i  i 


(20) 


and  the  subscripts  on  T  identify  one  of  the  six  orientations  considered. 
The  factor  3  in  front  of  the  expression  for  D  permits  D  to  reach  1.0 

XX  XX 

when  V  reaches  1.0  with  all  T’s  equal.  The  0.5  in  the  expression  for  D 

_  xx 

2  o 

is  cos  45  .  The  damage  quantity  in  the  y  and  z  directions  are  computed 
by  Eq.  19  with  a  natural  permutation  of  the  subscripts.  For  the  xy 
direction  the  damage  is 

D  =  3T  &y  [T  +  0.5(T  +  T  )]  .  (21) 

xy  F  xy  xx  yy 


T 


With  the  foregoing  definitions  of  danage  D,  we  can  examine  in 


more  detail  the  relaxation  process  in  Eq.  13.  That  equation  must  be 
modified  for  compression  where  we  can  expect  the  deviators  to  reduce 
with  damage  but  the  pressure  should  not  be  affected.  Toward  this  end, 
the  gross  pressure  P  is  related  to  the  solid  stresses  (X8^  for  the 
tensile  case.  Then  the  gross  stress  0^  is 


°ii  -  °Ii  (1  -  Dii> 


(22) 


P  = 


Saiia  -  V 


£cr'S(l  -  D  )  +  P  [3  -  (D  +  Dnn  +  D  )  ] 
ii  ii  s  11  22  33 


=  \  Ea,s  (1  -  D  )  +  P  (1  -  D) 
3  ii  ii  s 


(23) 


where  T,  D  =  3D,  P  is  the  solid  pressure,  O  is  stress  on  the 
x  x  s  i  A 

/  s 

gross  section  and  (T^  is  deviator  stress  on  the  intact  or  solid  material. 
For  the  compression  case  a  suitable  modification  of  Eq.  23  is  obtained 
by  deleting  the  factor  (1  -  D)  and  defining  a  modified  pressure  P7; 


»'  =  h  Ecr'fu  -  D  >  + 


P’  »  -  Eflr:.<i  -  d  >  +  p  =  -  Sa8  (l  -  d  )  +  p  d 
3  ii  ii  s3  ii  ii  s 


(24) 


=  P  +  P  D 
s 


Then  for  tension  let  P7  =  P  and  use  Eq.  24  to  define  P7  in  compression. 
With  these  definitions  of  P7,  the  solid  stresses  are  related  to  the 
gross  stresses  as  follows  in  both  tension  and  compression: 


_  /  _/ 
"ii  *p 


ii 


1  -  D 


(25) 


ii 
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At  the  beginning  of  a  strain  increment  calculation,  T  and  P  are 

s 

computed  from  the  existing  shear  bands  and  density,  and  P  is  calculated 

from  the  old  gross  stresses,  (X  and  from  P  for  compression.  Then  the 

H  8 

8 

or  are  calculated  from  Eq.  25.  At  the  end  ojf  the  increment,  Eqs.  24 
and  25  are  again  used  to  obtain  the  gross  stresses. 

The  criterion  for  full  fragmentation  is  given  in  Eq.  18.  This 
criterion  is  reached  simultaneously  with  the  condition  that  all  gross 
deviator  stresses  have  relaxed  to  zero. 

4.  Stress-Strain  Relations 

In  preparation  for  the  stress  calculations,  the  gross  stresses 
from  the  calling  program  are  transformed  to  stresses  in  the  intact 
material  using  Eq.  13.  These  stresses  are  then  separated  into  pressure 
and  deviator  components.  The  pressure  is  given  by  the  Mie-Gruneisen 
equation : 


P 


V1  '  5*’  *  rpE 


(26) 


where  p.  =  p/p  -  1 


r  is  Gruneisen  ratio 


p  is  the  current  density 


p  is  the  zero  pressure  density 


E  is  the  internal  energy 
2  3 

and  P  =  K  Up  +  K >  +  K  u.  is  the  pressure  on  the  Hugoniot.  K  is 
H  1  Z  3  1 

the  bulk  modulus,  and  K  and  K  are  coefficients  in  the  Hugoniot 

2  3 

expansion. 


The  deviator  stresses  are  computed  from  elastic,  plastic,  work¬ 
hardening,  rate-dependent  relations.  For  the  elastic  case  the  deviator 
stresses  are  computed  as 
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(27) 


/  .  /  / 
a  =  2GAc  +  (T 
ij  ij  ijo 


where  G  is  the  shear  modulus, 

is  the  deviator  strain  increment  and 

<3*  is  the  deviator  stress  on  the  previous  cycle, 
ijo 

For  yield  calculations,  the  effective  stress  from  von  Mises  plasticity 
theory  is  used 


a  =  J- £  a  a 
12  r,j  ij  i 


Yield  occurs  for  (J  2  Y,  the  yield  strength. 


A  linear  work-hardening  is  assumed  such  that 


Y  =  Y  +  Y  (P 
O  D 


where  Y  is  the  work -hardening  modulus 
D 


■/VffT 


P  P 

d<  dC  ,  the  equivalent  plastic  strain, 

i  j  i  j 


defined  in  terms  of  the  plastic  strain  deviator  increments  dC^  • 

Strain-rate  effects  are  treated  with  the  linear  viscous,  one- 


parameter  model 


d<J/  d€# 

-M  =  2G  -ii 

dt  dt 


dC 

-Y-ii 
iJ  dc-” 


where  T  is  a  time  constant.  This  equation  is  presumed  to  govern  only 
when  the  stresses  exceed  yield.  When  work-hardening  is  combined  with 
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the  rate  effects,  the  solution  must  be  obtained  numerically  for  general 

2 

cases.  The  solution  for  stress  derived  by  Seaman  et  al.  is 


a  =  O'  exp[-(t  -  t  )/t] 
ij  ijo  o 


(31) 


y5n 

Y  +  — —  0* 

1  4G  ij 


Y 


-N 

a 


— - —  (o1  N  -  a1  ) 

t  -  t  \  ij  ij©; 


1  -  exp[-(t  -  t  )/t] 
o 


where  t  -  t  is  the  duration  of  the  time  increment 
o 

Y  is  the  yield  strength  at  the  beginning  of  the  increment 
/  N  -N 

a  and  Or  are  elastically  computed  deviator  stress  and 

J 

effective  stress  from  Eqs.  27  and  28.  The  plastic  deviator  strain 
increments  govern  the  nucleation  and  growth  processes.  These  strains 
are  computed  as  the  difference  between  the  total  deviator  strains  and 
the  elastic  values. 

B.  Simulations  of  Fragmenting  Round  Experiments 

Computer  simulations  were  made  for  three  fragmenting  cylinders  of 
HF-1  high  fragmentation  steel,  one  cylinder  of  Armco  iron,  and  one  HF-1 
fragmenting  round  of  standard  munition  geometry.  The  shear  band  sub¬ 
routine  SHEAR2  was  used  as  the  stress-strain  relation  in  all  cases.  Four 
of  the  cases  were  considered  to  be  cylindrically  symmetric,  and  there¬ 
fore  calculations  were  made  with  the  one-dimensional  wave  propagation 
code  PUFF.  Since  the  fifth  round  (of  HF-1)  had  the  shape  of  a  standard 
munition,  a  two-dimensional  wave  propagation  code,  TROTT,  was  used  for 
that  case.  Adjustments  were  made  in  the  two  codes  so  that  the  identical 
SHEAR2  subroutine  was  used  with  each.  The  dimensionality  of  the  codes 
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refers  to  the  number  of  directions  along  which  quantities  can  vary 
(radial  in  PUFF,  radial  and  axial  in  TROTT) .  In  both  cases  the 
stresses,  strains,  and  damage  have  components  in  all  three  dimensions. 

The  geometries  of  the  five  experiments  are  shown  in  Figures  16  and 

4 

17.  Three  are  configurations  used  by  Crowe  et  al.;  the  other  two  are 
part  of  this  work.  Experiment  5084-3  was  simulated  first  with  TROTT 
and  SHEAR2  (with  no  damage  permitted),  and  the  pressure  history  at  mid¬ 
length  on  the  inner  radius  of  the  round  was  fitted  to  an  exponential. 
This  exponential  pressure  history  was  used  for  the  SRI  PUFF  simulations 
of  that  experiment.  For  all  other  calculations,  the  explosive  was 
simulated  as  a  detonating  material.  In  the  cylindrically  symmetric 
cases  (SRI  PUFF  calculations),  the  explosion  was  treated  as  a  simultane¬ 
ous  detonation.  For  the  TROTT  problem,  (the  fragmenting  projectile) 
the  explosive  was  simulated  with  a  running  detonation.  All  explosives 
were  treated  as  polytropic  gases  except  the  PETN  in  5084-4,  which  was 

handled  with  a  tabular  (pressure-volume)  equation  of  state.  The  tabular 

5 

values  were  obtained  from  a  TIGER  detonation  calculation.  More  details 
about  these  calculations,  including  the  data  decks  for  the  calculations, 
are  given  in  Appendix  B. 

The  results  of  the  simulation  of  5084-3  were  simply  that  the  maxi¬ 
mum  plastic  shear  strain  was  1.15  times  the  threshold  value.  Because 
the  experiment  showed  only  slight  evidence  of  the  beginning  of  shear 
banding,  there  seems  to  be  a  satisfactory  correlation  of  test  and 
calculation. 

The  fragment  size  distributions  for  the  three  cylinders  that  showed 

complete  separation  are  shown  in  Figure  18.  The  computed  results  are 

also  shown.  The  two  tests  of  Crowe  occurred  at  a  strain  rate  of  about 
4 

2  x  10  /sec  through  the  period  of  high  damage  while  experiment  5084-4 

3 

occurred  at  a  strain  rate  of  6  x  10  /sec.  The  experimental  results 
with  HF-1  show  that  an  increase  in  strain  rate  leads  to  a  much  larger 
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CUMULATIVE  NUMBER  OF  FRAGMENTS  GREATER  THAN  R,  PER  cm 


10 


MA-5084-28 


FIGURE  18  COMPARISON  OF  EXPERIMENTAL  AND  COMPUTED  FRAGMENT 
SIZE  DISTRIBUTIONS  FOR  HF-1  AND  ARMCO  IRON  CYLINDERS 


number  of  fragments,  and  consequently  to  smaller  fragments.  The  Armco 
iron  test,  while  conducted  at  the  high  strain  rate,  shows  a  fragment 
size  distribution  more  like  that  from  experiment  5084-4,  which  was 
conducted  at  the  lower  strain  rate.  This  difference  may  be  associated 
with  the  smaller  number  of  inclusions  in  Armco  iron,  or  with  the  lower 
yield  strength;  the  data  base  is  not  yet  adequate  to  provide  such  an 
association. 

The  marked  difference  in  the  distributions  measured  in  the  two 

HF-1  experiments  provide  a  good  test  of  the  model.  The  same  parameters 

were  used  in  both  HF-1  calculations.  The  correspondence  between  computed 

and  observed  distributions  Indicates  that  the  model  can  account  for 

the  differences  in  strain  rate,  geometrical  layout,  and  explosive 

pressure.  The  Armco  iron  calculation  was  run  with  a  nucleation  co- 

-3  -4 

efficient  one-third  that  for  HF-1  (C  =  10  for  HF-1,  3  x  10  for 

n 

Armco).  This  difference  in  parameters  probably  arises  because  of  the 
fewer  nucleation  sites  in  Armco  iron.  The  nucleation  plastic  strain 
threshold  probably  should  also  have  been  changed  for  Armco  iron,  but 
no  value  was  known. 

The  PUFF  simulations  of  the  three  cylinder  experiments  showed 
that  the  contained  rounds  required  much  longer  to  reach  full  fragmenta¬ 
tion  than  the  unconfined  rounds.  Crowe's  experiments  took  about  45  visec 

4 

to  reach  full  fragmentation,  according  to  the  calculations,  while  the 
contained  cylinder  required  over  200  paec.  In  all  three  cases  the 
shear  banding  began  at  the  inner  radius  and  gradually  moved  out.  All 
computational  cells  showed  shear  banding  in  all  three  of  the  active 
orientations  illustrated  in  Figure  8.  However,  in  each  case,  orienta¬ 
tion  1  was  the  first  one  activated;  the  others  began  to  nucleate  bands 
much  later.  When  the  bands  from  the  three  orientations  were  added, 
orientations  2  and  3  contributed  only  a  few  percent  at  most.  This 
result  agrees  with  the  experiments,  which  indicate  that  most  bands 
are  in  the  first  orientation. 
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The  two-dimensional  simulation  of  Crowe's  sample  fragmenting  pro¬ 
jectile  was  performed  with  the  cell  layout  shown  in  Figure  19.  All  the 
metal  parts  are  of  HF-1  steel.  No  provision  was  made  for  the  joint  in 
the  two  parts  of  the  case.  The  primary  explosive  in  the  fuse  is  CH6  and 
the  secondary  explosive  is  PBXN  106.  The  simulation  includes  treatment 
of  a  running  detonation  that  starts  in  the  primary  and  sweeps  at  the 
C-J  velocity  through  the  secondary  explosive.  The  progress  of  the  run¬ 
ning  detonation  is  evident  from  the  cell  motion  in  Figure  20,  which  shows 
a  series  of  plots  from  the  two-dimensional  calculations.  The  indicators 
in  the  figure  show  the  attendant  progress  of  fracturing  in  the  HF-1  cells. 

Figure  21  shows  a  comparison  of  the  computed  fragment  size  distri¬ 
bution  with  the  measured  distribution.  The  computed  distribution  was 
obtained  by  summing  the  fragment  size  distributions  from  all  computational 
cells,  accounting  for  the  relative  mass  associated  with  each  cell.  Cells 
that  did  not  fragment  were  mainly  in  the  nose,  and  all  unfragmented  masses 

were  added  in  as  a  single  large  chunk.  In  this  two  dimensional  simulation, 

o 

the  first  shear  band  mode  activated  was  usually  orientation  1  at  45 
between  radial  and  circumferential  directions.  Usually  this  mode  also 
dominated  at  the  time  of  fragmentation.  Considerable  shear  strain  also 
occurred  along  planes  of  constant  radius  or  axial  distance,  and  if  bands 
had  been  permitted  on  these  orientations,  some  shear  bands  would  have 
occurred  in  the  calculations.  The  duration  of  detonation  was  about  65 
psec.  By  110  psec,  nearly  all  the  shear  banding  cells  were  completely 
fragmented.  No  shear  banding  occurred  forward  of  K-row  48. 

Having  completed  these  simulations  of  the  fragmenting  rounds,  we 
then  examined  ancillary  aspects  of  the  computed  results  to  see  how 
these  aspects  matched  the  data.  Among  these  aspects  were  the  amount  of 
homogeneous  plastic  strain,  the  effect  of  microcracking,  the  relative 
importance  of  the  nucleation  process,  and  the  information  that  can  be 
obtained  from  the  shape  of  the  fragment  size  distribution. 
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•  Some  Shear  Banding 
■  Fragmentation  by  Shear  Banding 


FIGURE  20  PROGRESS  OF  DETONATION  AND  SHEAR  BAND  DAMAGE 
COMPUTED  IN  THE  TWO-DIMENSIONAL  SIMULATION 
OF  CROWE'S  SAMPLE  FRAGMENTING  PROJECTILE 
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MA-5084-40 

FIGURE  21  COMPARISON  OF  COMPUTED  AND  OBSERVED  FRAGMENT  SIZE 

DISTRIBUTIONS  FOR  CROWE'S  SAMPLE  FRAGMENTING  PROJECTILE 
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In  the  model  calculation*,  It  was  assumed  that  the  shear  band 
nucleation  and  growth  gradually  takes  a  larger  portion  of  the  total 
plastic  strain  until  all  is  being  absorbed  by  the  bands.  This  means 
that  after  the  threshold  of  30%  equivalent  plastic  strain  is  reached, 
the  plastic  strain  taken  in  homogeneous  deformation  will  continue 
to  increase  until  fragmentation  occurs.  In  our  computed  results  the 
total  equivalent  plastic  strain  imposed  is  60  to  80%.  But  the  homogeneous 
plastic  strain  at  the  time  of  fragmentation  is  about  40% — the  rest  is 
used  in  producing  bands.  Experimental  measurements  indicate  that  there 
was  only  35%  homogeneous  strain  in  the  fragments.  It  may  be  necessary 
to  reexamine  the  nucleation  and  growth  processes  to  find  means  whereby 
the  band  formation  mechanism  can  absorb  a  larger  fraction  of  the  induced 
plastic  strain. 

In  the  calculations,  fragmentation  was  assumed  to  occur  only  by 
shear  banding,  whereas  in  some  of  the  experiments  microcracking  in  the 
outer  layer  participated  in  the  fragmentation  process.  In  future  work 
the  effects  of  microcracking  must  be  included  in  the  computations. 

The  computed  fragment  size  distributions  were  obtained  by  varying 
mainly  the  nucleation  rate  coefficient.  The  Orowan  equation  appears 
to  force  a  balance  between  nucleation  and  growth  so  that  an  increase  in 
nucleation  automatically  results  in  less  growth. 

The  fragment  size  distributions  generally  are  linear  on  a  semilog 
graph  of  cumulative  number  versus  fragment  radius.  However,  there  are 
some  interesting  deviations  from  linearity  in  the  calculations  and  the 
experiments.  A  concave  upward  portion  for  small  radii  results  from 
the  nucleation  of  more  bands  at  a  late  stage.  A  concave  downward 
portion  at  the  large  radius  end  results  from  reaching  a  fixed  velocity 
limit  such  as  the  shear  wave  velocity. 
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C.  Conclusions  and  Recommendations 


A  constitutive  model  was  developed  to  provide  elastic-plastic 
stress-strain  relations  and  to  represent  the  shear  banding  phenomena 
observed  in  a  series  of  fragmentation  experiments  in  HF-1  steel.  The 
model  includes  the  nucleation  of  shear  bands,  growth  of  bands,  and 
coalescence  of  bands  to  form  fragments.  The  model  accounts  for  stress, 
strain,  and  damage  quantities  in  three  dimensions,  although  it  is  used 
in  one-  and  two-dimensional  wave  propagation  codes.  The  model  permits 
bands  to  grow  in  six  specific  orientations;  these  include  the  three 
that  can  occur  in  exploding  cylinders  plus  three  more  that  may  be 
activated  in  impacts.  The  developing  damage  causes  relaxation  of  the 
stress  in  directions  normal  to  the  planes  of  the  bands;  hence,  growth 
of  bands  develops  an  anisotropy  in  the  material. 

The  results  of  calculations  with  the  shear  band  model  include  stress 
histories,  and  the  size  and  number  of  bands  in  each  orientation.  The 
band  distributions  can  be  transformed  to  a  fragment  size  distribution 
and  summed  over  all  the  material  to  form  a  Mott  plot. 

Five  experiments  were  simulated:  four  of  HF-1  steel  and  one  of 
Armco  iron.  Three  cylinders  of  HF-1  showed  good  correspondence  between 
experimental  and  computed  fragment  size  distributions  although  the 
strain  rates,  geometry,  and  explosive  power  were  significantly  different 
for  the  three  cases.  The  shell  in  the  shape  of  a  standard  round  of  HF-1 
and  the  Armco  iron  cylinder  also  showed  good  correspondence  between 
experimental  and  computed  fragment  size  distributions.  The  model  gave 
the  correct  orientation  for  most  of  the  bands,  and  partitioned  the 
plastic  strain  between  homogeneous  strain  and  that  absorbed  during 
shear  banding.  These  results  Indicate  that  the  nucleation  and  growth 
processes  of  the  model  are  appropriate  and  are  able  to  account  for 
strain  rate,  geometry,  and  pressure  differences. 
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The  present  shear  band  model  was  developed  for  the  current  experi¬ 
mental  data  on  HF-1  steel,  and  has  been  verified  only  for  those  data. 
These  data  are  for  cylindrical  geometry  which  is  not  sensitive  to  many 
features  of  the  model:  hence,  those  features  are  not  constrained  by 
the  data.  To  develop  a  model  representing  a  wide  range  of  material 
types,  loading  conditions,  and  geometries  will  require  more  experiments 
and  coordinated  model  development.  Specific  areas  of  the  model  that 
are  recommended  for  further  study  are  the  following: 

(1)  The  nucleation  law  requires  verification  experimentally,  or 
modifications  to  match  experiments.  Associated  experiments 
should  be  conducted  for  a  broad  range  of  strain  rates  and 

for  various  levels  of  damage.  Especially  important  are  experi¬ 
ments  at  the  level  of  incipient  damage,  because  at  that  point 
the  nucleation  process  has  been  little  affected  by  the  develop¬ 
ment  of  damage. 

(2)  The  growth  process  should  be  examined  in  the  light  of  the 
experimental  data  obtained  in  (1)  above.  The  low-damage 
experiments  should  guide  selection  of  a  new  function,  if 
necessary,  and  evaluation  of  the  parameters.  Band  motion  may 
also  be  examined  analytically  by  studying  the  behavior  of  a 
single  band  in  homogeneous  material.  Band  width  should  be 
determined  as  a  function  of  strain  rate,  thermal  conductivity, 
thermal  strength  reduction,  and  work  hardening.  These  width 
calculations  can  probably  be  performed  with  a  one-dimensional 
shear  wave  propagation  code.  Growth  of  the  bands  may  be 
treated  in  the  same  way  with  a  two-dimensional  code.  If  both 
width  and  growth  velocity  are  governed  by  homogeneous  processes, 
these  computational  studies  could  aid  greatly  in  defining  some 
of  the  micro-mechanisms  of  the  shear  band  process. 
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(3)  In  Impacts  and  fragmenting  rounds,  brittle  fracture  and  shear 
bands  combine  to  form  fragments,  and  this  combined  process  is 


not  yet  incorporated  into  the  NAG  (Nucleation  and  Growth) 

* 

models.  For  this  process,  fragments  as  well  as  material  at 
incipient  damage  levels  should  be  examined.  Then  analytical 
techniques  for  joining  the  two  damage  processes  should  be 
developed . 

(4)  In  the  model,  both  nucleation  and  growth  in  any  direction  are 
presumed  proportional  to  the  plastic  shear  strain  in  the  rele¬ 
vant  direction.  This  hypothesized  relationship  should  be 
examined  in  the  light  of  data  from  incipient  damage  experiments. 


The  combined  process  was  accounted  for  in  the  work  reported  in  Ref.  2 
by  permitting  cells  to  undergo  either  brittle  fracture  or  shear  banding, 
but  not  both.  In  that  study  only  one  orientation  of  shear  bands  was 
permitted. 
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V  OBLIQUE  ROD  IMPACT  -  EXPERIMENTAL  PHASE 


The  mechanisms  for  fracture  of  a  long-rod  penetrator  under  high 
velocity  impact  (either  oblique  or  perpendicular)  with  an  armor  plate 
are  inadequately  understood.  Large  pressures  and  large  shear  deformations 
at  high  strain  rates  make  this  a  difficult  problem  to  analyze  or  model. 

In  a  preliminary  effort  to  establish  the  key  penetrator  failure 
modes,  we  performed  the  following  tasks: 

(1)  Examined  metallographically  the  distal  ends  of  steel  and 
tungsten  alloy  long-rod  penetrators  recovered  from 
experiments  at  Lawrence  Livermore  Laboratories  (LLL) . 

(2)  Examined  metallographically  fragments  that  LLL  recovered 
from  the  proximal  end  of  obliquely-impacted  tungsten  alloy 
rods . 

(3)  Performed  a  series  of  oblique-angle  long-rod  penetrator 
experiments  using  HF-1  steel  to  study  break-up  characteristics 
as  a  function  of  velocity. 

(4)  Examined  the  fragments  recovered  from  a  long-rod  tungsten 
alloy  penetrator  perpendicular  impact  performed  at  the 
Ballistic  Research  Laboratory. 

These  activities  are  described  in  detail  in  the  following  pages. 

A.  LLL  Experiments 

From  Dr.  John  Scudder  and  Mr.  Leroy  Hoard  of  LLL,  we  obtained 

fragments  recovered  from  scale-model  long-rod  penetrator  oblique-impact 

tests  performed  in  their  laboratories.  The  projectiles  all  had  a 

diameter  of  3.9  mm,  a  length  of  39  mm,  and  a  hemispherical  nose,  and 

o 

were  impacted  at  a  velocity  of  «  1  km/sec  and  a  55  angle  of  obliquity 

* 

into  a  target  of  mild  steel  with  an  overlay  of  AD85  alumina  ceramic. 


Product  of  Coors  Porcelain  Co.,  Golden,  Colorado. 
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We  first  examined  the  distal  halves  of  the  projectiles - the  only 

fragments  large  enough  to  be  recovered  without  any  special  recovery  tech¬ 
niques.  Figure  22(a)  shows  a  metallographic  specimen  made  by  cutting 
through  the  axis  of  the  specimen  parallel  to  the  plane  of  greatest 
bending  (all  subsequent  specimens  were  cut  in  the  same  way).  The  material 
is  Saratoga  AISI-01  oil  hardened  tool  steel  (Rc  48) .  Photomicrographs 
of  a  region  near  the  forward  edge  of  the  specimen  show  clearly  a  shear 
band  intersecting  the  edge  and  a  tensile  crack  extending  part  way  down  the 
length  of  the  shear  band.  Since  this  was  the  only  internal  crack  visible 
in  this  specimen  and  since  there  was  also  a  transformed  region  (indicating 
a  shear  band)  along  the  fractured  forward  edge  of  the  specimen,  there  is 
evidence  that  shear  banding  and  subsequent  tensile  failure  along  a  shear 
band  is  an  important  failure  mechanism  for  a  steel  long-rod  penetrator, 
at  least  in  the  distal  regions. 

We  next  examined  the  distal  halves  of  tungsten  alloy  projectiles. 

The  alloy  contained  90%  tungsten,  7%  nickel,  and  3%  iron  by  weight,  and 
was  composed  of  nearly  spherical  tungsten  grains  a-  5  u,m  in  diameter  in 
a  matrix.  Unlike  steel,  tungsten  does  not  exhibit  a  readily  observed 
phase  change  under  high  temperatures,  which  can  be  used  as  a  clear 
indicator  of  shear  banding.  However,  this  tungsten  alloy  does  have  an 
advantage  over  some  steel  in  that  strain  deformation  is  indicated  by  grain 
elongation.  Figure  23  shows  three  photomicrographs  from  a  tungsten 

penetrator - one  taken  from  an  undeformed  region,  showing  largely  circular 

grain  cross  sections;  one  from  a  region  of  large  but  primarily  homogeneous 
shear  deformation,  showing  grain  cross  sections  with  aspect  ratios  of 
around  4  to  1;  and  one  from  a  region  along  the  proximal  edge  where 
inhomogeneous  shear  deformations  have  elongated  grains  to  aspect  ratios 
as  high  as  20  to  1,  and  caused  shear  cracking. 
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FIGURE  22  SPECIMEN  (a)  FROM  DISTAL  END  OF  OBLIQUELY  IMPACTED  TOOL  STEEL 
LONG-ROD  PENETRATOR  FIRED  AT  LAWRENCE  LIVERMORE  LABORATORY 
WITH  PHOTOMICROGRAPHS  (b)  AND  (c)  SHOWING  SHEAR  BANDING 


MP-5084-19 


FIGURE  23  PHOTOMICROGRAPHS  OF  DISTAL  HALF  OF  OBLIQUELY  IMPACTED  TUNGSTEN 
ALLOY  LONG-ROD  PENFTRATOR  SHOWING  (al  UNDEFORMED  REGION. 

(b)  REGION  OF  PRIMARILY  HOMOGENEOUS  SHEAR  DEFORMATION,  AND 

(c)  REGION  OF  LARGE  INHOMOGENEOUS  SHEAR  DEFORMATION  RESULTING 
IN  SHFAR  CRACK 


We  have  thus  seen  evidence  of  failure  due  to  shear  banding  or 
inhomogeneous  shear  deformations  in  the  distal  region  of  both  steel 
and  tungsten  alloy  long-rod  penetrators.  Since  the  shear  strains  and 
strain  rates  would  be  higher  near  the  proximal  ends,  this  type  of  failure 
mechanism  should  be  even  more  important  there.  To  test  this  belief,  we 
suggested  to  Dr.  Scudder  and  Mr.  Hoard  a  means  of  recovering  fragments 
from  the  proximal  halves  of  obliquely  impacted  long-rod  penetrators  by 
capturing  them  in  a  relatively  soft  rubber-like  substance.  They  subsequently 
performed  such  experiments,  and  we  examined  several  recovered  tungsten  alloy 
specimens  one  of  which  is  shown  in  Figure  24.  The  structure  is  very  similar 
to  that  shown  in  Figure  23.  In  Figure  24(b),  a  gradual  transition  from  no 
shear  deformations  to  large  shear  deformations  takes  place  across  the  width 
of  the  specimen.  And  in  Figure  24(c),  narrow  bands  of  extremely  high  shear 
deformation  are  seen  between  regions  of  moderate  shear  deformation.  We 
cannot  know  precisely  where  in  the  penetrator  this  fragment  originated, 
but  it  was  within  the  proximal  half. 

Our  conclusion  is  that  failure  by  shear  banding  or  inhomogeneous  shear 
deformation  occurs  over  a  large  region  in  a  long-rod  penetrator. 

B.  SRI  Experiments 

SRI's  smooth-bore-gun  facility,  which  had  been  used  in  previous  armor 

penetration  experiments,^  was  modified  for  use  in  scale-model  long-rod 

penetrator  oblique  impact  experiments.  Hemispherical-nosed  HF-1  steel 

o 

projectiles  3.175  mm  in  diameter  and  38.1  mm  long  were  fired  at  a  55 
angle  of  obliquity  into  1-cm-thick  steel  armor  plates.  The  relatively 
thick  armor  was  chosen  so  as  to  minimize  damage  to  the  target,  and  thereby 
make  the  projectile  failure  easier  to  model  computationally.  HF-1  steel 
was  chosen  because  it  is  a  material  for  which  dynamic  fracture  parameters 
are  known. 
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(a) 


(b) 


(0 

MP-5084-17 


FIGURE  24  FRAGMENT  (a)  RECOVERED  FROM  OBLIQUELY  IMPACTED  TUNGSTEN  ALLOY 
LONG-ROD  PENETRATOR  FIRED  AT  LAWRENCE  LIVERMORE  LABORATORY 
WITH  PHOTOMICROGRAPHS  SHOWING  EVIDENCE  OF  LARGE  HOMOGENEOUS 
(b)  AND  NONHOMOGENEOUS  (c)  SHEAR  DEFORMATIONS 
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We  fired  a  series  of  shots  at  impact  velocities  ranging  from  0.17 
to  0.97  km/sec,  and  the  results  are  shown  in  Figure  25.  As  the  impact 
velocity  increases  to  about  2/3  km/sec,  the  projectiles  become  more  and 
more  bent,  eventually  separating  into  several  pieces.  However,  at  higher 
velocities,  the  bending  decreases  and  the  projectiles  undergo  large  cross- 
sectional  deformations  in  their  proximal  regions.  In  all  cases,  target 
damage  was  negligible. 

The  damaged  projectiles  were  sliced  in  half  along  their  plane  of 
maximum  bending,  polished,  etched,  and  examined  metallographically . 

The  surfaces  of  separation  appeared  to  result  from  brittle,  or  tensile 
failure,  but  there  was  no  evidence  of  either  shear  banding  or  tensile 
failure  in  the  proximal  region.  The  higher  velocity  shots  exhibited 
a  transformed  region  along  the  outer,  flattened  surface,  which  undoubtedly 
resulted  from  the  projectile  sliding  along  the  surface  of  the  target  as 
it  was  being  deflected. 

C.  BRL  Experiments 

Late  in  the  program  Dr.  Zukas  of  Ballistics  Research  Laboratory 

sent  us  material  recovered  from  a  long-rod  penetrator  impact  experiment 

performed  at  BRL.  The  material  was  the  same  tungsten  alloy  used  in  the 

LLL  penetrators;  however,  the  scale  was  larger — the  hemispherical-nosed 

projectile  was  0.79  cm  in  diameter  and  7.9  cm  in  length,  and  the  impact 

o 

occurred  at  0.850  km/sec  with  a  0  angle  of  obliquity  into  a  2.54-cm-thick 
steel  target. 

The  recovered  fragments  included  the  distal  quarter  and  an  assortment 
of  small  fragments.  Figure  26  shows  a  slice  through  the  large  fragment, 
with  a  close-up  of  the  area  that  contains  what  appears  to  be  a  long  partially 
separated  shear  band.  Detailed  photomicrographs,  seen  in  Figure  27,  prove 
that  it  is  indeed  a  shear  band.  The  high  magnification  insets  clearly 
show  bands  less  than  10  M>m  wide,  across  which  the  two  opposite  surfaces 
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55°  IMPACT  HF-1  PROJECTILES 


0.97  mm/jisec 


0.84  mm/Msec 

0.61  mm//i$ec 
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0.30  mm/jusec 


0.17  mm/psec 
VIRGIN 


MP-5064-18 

FIGURE  25  HF-1  STEEL  LONG-ROD  PROJECTILES 
RECOVERED  FROM  OBLIQUE  IMPACT 
TESTS  AT  SRI  INTERNATIONAL 
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FIGURE  26  POLISHED  SECTION  (i 

LONG-ROD  PENETRA' 
BAND 
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have  slid  as  much  as  50  pm  or  more.  The  inset  at  lower  left  shows  two 
such  bands  running  parallel  to  each  other  about  60  pm  apart. 

A  polished  section  through  one  of  the  smaller  fragments  is  seen  in 
Figure  28,  along  with  two  high  magnification  insets  revealing  shear- 
banded  regions.  These  particular  bands  are  quite  curved,  which  indicates 
either  that  the  material  was  in  a  highly  nonuniform  state  of  shear  strain 
prior  to  the  banding  or  that  the  material  continued  to  undergo  large 
homogeneous  shear  deformations  after  the  band  formed.  Also,  the  highly 
deformed  grain  that  exists  away  from  the  banded  regions  is  indicative  of 
the  toughness  of  this  alloy,  which  will  sustain  large  amounts  of  shear 
strain  before  banding. 

D.  Conclusions 

Shear  banding  is  a  principal  mode  of  material  failure  in  high-velocity 
Impacts  of  tungsten  alloy  and  steel  long-rod  penetrators,  although  the  SRI 
experiments  Indicate  that  tensile  fracture  may  also  be  important,  particularly 
at  larger  impact  angles. 
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VI  OBLIQUE  ROD  IMPACT-COMPUTATIONAL  PHASE 


Oblique  impact  calculations  were  conducted  with  a  two-dimensional 

wave  propagation  computer  program.  The  calculations  simulate  a  plate 

impacting  on  edge  another  plate  at  an  oblique  angle.  Both  the  conditions 

in  which  the  projectile  sticks  to  the  target  and  in  which  it  slips  were 

o  o  o 

treated.  For  each  condition.  Impact  angles  of  30  ,  45  ,  and  55  to 
normal  were  considered. 

The  purpose  of  the  calculations  was  to  examine  the  phenomena 
associated  with  oblique  impacts  to  gain  insight  into  the  reasons  for 
breakup  of  the  projectile  and  the  mechanism  of  damage  in  the  impact 
region.  We  examined  the  regions  of  fracture  damage  and  locations  near 
the  point  of  impact  where  shear,  compression,  and  bending  stresses 
occurred,  and  we  noted  how  these  locations  changed  with  the  degree  of 
slip  along  the  target.  We  used  an  active  fracture  damage  model  (with 
stress-relaxation  associated  with  damage)  to  see  whether  the  damage 
would  significantly  affect  subsequent  bending  and  projectile  breakup. 

The  computer  program  used  is  called  TROTT.  It  is  a  two-dimensional, 
finite-difference  Lagrangian  code,  and  is  similar  to  HEMP  and  TOODY. 

Three  features  made  it  attractive  for  this  problem: 

(1)  TROTT  is  small;  all  storage  is  in  core. 

(2)  The  array  storage  is  arranged  to  permit  very  large  data 
storage  for  cells  undergoing  fracture. 

(3)  The  code  is  simple  and  easy  to  modify  for  the  special 
impact  boundary  conditions. 

During  the  project,  modifications  were  made  to  allow  improved  treatment 
of  the  impact  problem. 


The  remainder  of  this  section  treats  the  special  modifications  made 

to  the  code  TROTT,  the  geometries  used  in  the  calculations,  and  the  results 
of  the  simulations. 

A.  Code  Modifications 

The  TROTT  code  was  altered  to  include  a  special  nonnormal  impact 
condition  and  triangular  cells.  The  projectile  is  divided  into  the  usual 
discrete  finite-difference  cells,  whereas  the  target  is  treated  as  rigid 
and  immovable.  The  rigid  target  was  used  for  simplicity  because  the  main 
interest  was  in  breakup  of  the  projectile  at  some  distance  back  from  the 
impact  plane.  It  was  felt  that  the  conditions  of  sticking  and  slipping 
would  bracket  the  boundary  conditions  actually  experienced  by  a  projectile. 

For  the  slipping  boundary  condition,  coordinate  motion  was  permitted 
as  usual  until  the  point  reached  the  target.  Thereafter  only  tangential 
motion  along  the  target  face  was  allowed.  For  example,  in  Figure  29, 

Point  A  moves  until  it  strikes  the  target  at  B  and  then  slips  to  point  D. 
Let  (x^,  y^)  represent  one  point  on  the  target.  Then  the  equation  of  the 
target  is 

x  =  x^  +  (y  -  y^)  tan  0  (33) 

where  0  is  the  angle  of  the  target  measured  so  that  0  =  0  is  a  normal 
impact.  Then  if  a  coordinate  with  the  initial  position  (x  ,  y  ),  that  is 

A  A 

point  A  in  Figure  29,  moves  with  velocities  x,  y,  impact  occurs  at  point 
B  at  a  time 

*  -  *  +  (y  -  y,>  tan  0 

fit  =  — - - - - - - -  (34) 

A  -  £  tan  0 

If  fit  is  less  than  the  current  time  step  in  the  wave  propagation 

calculations  (i.e.,  point  C  is  beyond  point  B) ,  then  the  impact 
* 

coordinate  locations  are  at  B 
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FIGURE  29  SLIDING  OR  STICKING  CONDITIONS  AT  OBLIQUE  TARGET 
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*B  ■  *A  ♦  4  8t 


(35) 


yB  "  yA  *  }  5t 

After  the  impact,  cell  motion  is  tangential.  The  tangential  velocity  is 

VT  =  x  sin  ft  +  7  cos  0  . 

The  x  and  y  coordinates  of  this  velocity  are 

V  =  V  sin  ft 

x  T 

V  =  T  cos  0. 

y  T 

The  final  coordinate  locations  are  then 

*D  =  *B  +  \  (At  -  6t) 

yn  *  y„  +  v  (At  -  6t) 

D  B  y 

where  At  is  the  current  time  step. 

The  sticking  boundary  condition  required  that  no  motion  occurs  after 
impact.  First,  the  coordinate  motions  were  computed  in  the  usual  way. 

If  a  point  had  penetrated  the  target  plane  (i.e.,  point  C  in  Figure  29), 
it  was  repositioned  at  the  point  where  it  first  touched  the  target, 
point  B,  using  Equation  35.  Velocities  of  the  impacting  coordinates 
were  set  to  zero. 

A  triangular  cell  feature  was  added  to  TROTT  to  handle  the  large 
distortion  problem  near  the  plane  of  impact.  The  usual  quadratic  cells 
become  excessively  distorted  in  the  large  shear  flows  near  the  front  and 
they  either  halt  the  calculation  or  require  rezoning  of  the  cells. 
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The  triangular  cells  are  very  resistant  to  cell  distortion  and  therefore 

can  follow  the  impact  somewhat  better.  The  triangular  cells  were  constructed 

by  simply  dividing  quadrilateral  cells  in  two,  and  computing  strain  and 

stress  separately  for  each  half.  The  analysis  followed  the  lines  suggested 
6 

by  Johnson  for  impact  calculations. 

Following  incorporation  of  the  triangular  cells,  we  performed  several 
short  calculations  with  only  a  few  cells  to  examine  the  new  properties. 

The  study,  though  not  exhaustive,  was  adequate  to  guide  our  use  of  the 
triangular  cells  and  to  interpret  the  results.  It  revealed  the  following: 

e  The  triangle  Q  artificial  viscosity  stresses  used 

with  quadrilateral  cells  are  also  needed  with  triangular 
cells  to  reduce  oscillations. 

e  The  time  steps  required  for  stability  of  triangular  cell 
calculations  are  the  same  as  those  of  the  quadrilateral 
cell: 

AX 

At  =  0.9  ~ 

C 

where  AX  is  the  length  of  the  shortest  side,  C  is  the 
sound  speed  augmented  to  account  for  artificial  viscosity, 
and  0.9  is  the  safety  factor. 

•  In  pure  compression  a  pair  of  triangular  cells  each  give 
stresses  identical  to  the  quadrilateral  cell  from  which 
they  were  made. 

•  In  shear,  one  of  the  triangular  cells  in  a  pair  goes  into 
tension,  the  other  into  compression,  with  the  average  of  the 
two  stresses  being  about  that  of  a  single  quadrilateral  cell. 

These  stresses  are  often  quite  large,  thus  giving  a  misleading 
picture  of  the  stress  state. 
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The  last  item  above  shows  that  only  simple  equations  of  state  should  be 
used  for  the  triangular  cells.  More  complex  constitutive  relations  with 
threshold  levels  and  extreme  nonlinearities  would  only  give  misleading 
results.  Hence,  in  our  impact  calculations  we  used  the  triangular  cells 
only  for  the  first  few  rows  of  cells  at  the  head  of  the  projectile  and 
permitted  only  elastic-plastic  behavior  in  them. 

B.  Simulations 

The  geometry  used  in  the  calculations  is  shown  in  Figure  30.  Two- 
dimensional  plane  strain  conditions  are  used.  The  projectile  slab  has  a 
rounded  nose  as  shown  so  that  one  side  of  one  cell  is  in  continuous  contact 
with  the  target  at  impact.  A  projectile  velocity  of  1  mm/p,sec  was  used 
in  all  cases.  The  projectile  slab,  which  is  Armco  iron,  is  treated  as  an 
elastic-plastic  material  in  compression  and  with  brittle  fracture  (BFRACT 
subroutine)  in  tension.  In  addition  to  the  usual  stress  and  motion  results, 
moments  and  forces  in  the  impacting  slab  were  computed. 

The  resulting  shapes  of  some  of  the  projectiles  after  14  jisec  is 
shown  in  Figures  31  to  36. 

Following  Impact,  the  projectile  is  gradually  deformed  so  that  it  is 
in  contact  with  the  target  over  the  entire  projectile  nose  and  for  some 
distance  along  the  side  nearest  the  target.  Projectile  deformations 
appear  as  high  shears  that  distort  the  impact-plane  cells  and  elongate 
the  cells  just  behind  the  impact.  The  cells  are  elongated  through  the 
thickness,  suggesting  a  tendency  toward  splitting.  Cells  on  the  side  of 
the  projectile  farthest  from  the  target  are  also  elongated  in  the  axial 
direction,  indicating  a  tendency  to  break  in  bending. 

In  all  cases,  the  projectile  tends  to  enlarge  radially.  This  expansion 
is  more  evident  for  the  shallower  angles  of  impact.  Associated  with  this 
is  fracture  in  the  radial  direction  by  cells  near  the  center  of  the  projectile 
and  axial  fracture  on  the  projectile  surface  farthest  from  impact. 
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FIGURE  30  GEOMETRY  FOR  SIMULATIONS  OF  PROJECTILE  PLATE  IMPACTING  A  WALL 
AT  AN  OBLIQUE  ANGLE 


cm 


MA-5084-32 


FIGURE  31  DISTORTION  AND  DAMAGE  NEAR  FRONT  OF  PROJECTILE  PLATE  FOR  30° 
OBLIQUITY,  NO  SLIPPING 
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FIGURE  33  DISTORTION  AND  DAMAGE  NEAR  FRONT  OF  PROJECTILE  PLATE  FOR  45‘ 
OBLIQUITY,  NO  SLIPPING 


cm 


FIGURE  34  DISTORTION  AND  DAMAGE  NEAR  FRONT  OF  PROJECTILE  PLATE  FOR  45° 
OBLIQUITY,  SLIPPING 
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•  FRACTURE  HAS  BEGUN 


MA-5084-36 


FIGURE  35  DISTORTION  AND  DAMAGE  NEAR  FRONT  OF  PROJECTILE  PLATE  FOR  55° 
OBLIQUITY,  NO  SLIPPING 
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FIGURE  36  DISTORTION  AND  DAMAGE  NEAR  FRONT  OF 
FIGURE  36  °8l|QU)TY  supp(NG 


Graphs  of  the  forces  and  bending  moments  on  sections  through  the 
projectile  at  14  y,sec  are  shown  In  Figures  37  and  38.  These  values  are 
computed  by  summing  stresses  from  cells  with  the  same  initial  Lagranglan 
position.  The  force  and  moment  distributions  build  up  gradually  with 
time  toward  the  configurations  shown  in  Figures  37  and  38,  but  they  are 
fairly  constant  by  14  p,sec.  A  comparison  of  forces  in  these  figures  shows 
that  the  forces  have  a  higher  amplitude  and  extend  deeper  into  the 

projectile  for  the  sticking  case.  The  moments  show  a  similar  trend  for  the 

o  o  o 

45  and  55  impacts,  but  a  contrasting  trend  is  shown  at  30  .  These 

moment  variations  should  be  studied  further  with  calculations  performed 

for  longer  durations.  The  results  also  shed  light  on  the  SRI  oblique 

impact  experiments  reported  in  the  previous  section.  Comparison  of 

Figures  32,  34,  and  36  show  that  the  shear  strain  in  the  projectile  nose 

o 

is  dramatically  less  for  the  55  (slipping)  impact  than  for  impacts  at 
lower  angles.  Thus,  shear  banding  should  become  less  prevalent  at  higher 
impact  angles. 

Future  work  should  include  the  shear  banding  model  as  well  as  the 
tensile  fracture  model,  and  would  help  us  gain  an  understanding  of  the 
interaction  of  these  competitive  failure  modes. 
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FIGURE  38  BENDING  MOMENTS  IN  THE  IMPACTING  PLATE  AS  A  FUNCTION  OF  DISTANCE  ALONG 
THE  PLATE 


APPENDIX  A 


EQUIVALENT  PLASTIC  STRAIN  FOR  A  RADIALLY  EXPANDING  CYLINDRICAL  TUBE 


Consider  a  cylindrical  tube,  loaded  by  an  internal  pressure, 
expanding  radially  outward  with  no  motion  in  the  axial  direction.  What 
is  the  equivalent  plastic  strain,  averaged  across  the  wall  thickness, 
that  the  cylinder  has  undergone  after  expansion  and  after  the  internal 
pressure  has  been  relieved? 

Consider  the  portion  of  a  cylindrical  tube  shown  in  Figure  39. 

Let  r  =  radius  of  the  inner  surface,  R  =  radius  of  outer  surface, 
a  =  R-r  =  wall  thickness,  and  x  =  £(R+r)=  radius  of  the  midpoint. 
Furthermore  let  subscripts  o  and  f  refer  to  the  original  and  final 
dimensions,  respectively. 

Our  assumptions  are  that  there  is  no  axial  strain  (ez  =  0)  and 

that  the  final  volume  strain  is  zero  Gg+  *  o) .  Therefore, 

e  =  ~ea,  and  furthermore,  the  cross-sectional  area  of  the  tube  remains 
r  o  . 

/  a,  x 

constant  /i.e.,  a  x  =  ax, ,  or  f  =  o 
I  o  o  f  f  —  — 


Now  we  define  the  equivalent  plastic  strain : 

*vr- 

where  the  superscript  P  on  the  strain  terms  denotes  plastic  strains, 

*  This  definition  of  equivalent  plastic  strain  coincides  with  the 

definition  of  the  differential  dfP  in  Eq.  4  if  the  ratios  deP/de 
P  .  P  0  x 

and  dC^/dC^  are  constant  throughout  the  flow.  For  C  =0,  the 

ratios  are  constant.  z 
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-  Final  Position 


MA-5084-22 

FIGURE  39  GEOMETRY  AND  NOMENCLATURE  FOR  EQUIVALENT 
PLASTIC  STRAIN  DERIVATION 
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which  are  equal  to  the  final  strains,  since  there  are  no  residual 
stresses.  The  strains  may  now  be  computed  based  on  the  change  of 
wall  thickness,  or  on  the  change  of  radius.  A  change  in  thickness 
causes  a  change  in  radial  strain,  der  *  da/a.  Then 


In 


*f  -  rf 
Ro  -  ro. 


Similarly  a  change  in  radius  alters  the  circumference  and  induces  a 


circumferential  strain  de 


dr/r . 


In  these  expressions  and  Eg  represent  average  strains  through  the  wall 
thickness.  The  equivalent  plastic  strain  is  then 


ep  -  /3|er| 
“  ^3 1  eQ  | 


For  intact  material  the  strains  based  on  thickness  or  radius  changes  are 
equivalent.  But  for  fracturing  material,  the  thickness  measurement 
gives  the  homogeneous  plastic  strain,  whereas  the  radius  gives  the  total 
plastic  strain  (homogeneous  plus  that  takes  by  the  bands). 
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APPENDIX  B 


USE  OF  THE  SHEAR  BANDING  SUBROUTINE  IN  WAVE 
PROPAGATION  CALCULATIONS 

The  shear  band  subroutine  SHEAR2  is  used  in  one-  and  two-dimensional 
wave  propagation  codes  to  predict  shear  banding  and  fracture  under  high 
rate  loadings.  The  procedure  for  inserting  the  subroutine  into  wave 
propagation  codes  is  described  below.  A  brief  outline  of  the  subroutine 
is  given,  but  the  basic  equations  are  derived  in  Chapter  IV.  The  shear 
banding  parameters  are  described  and  some  guidance  is  given  for  choice 
of  values.  The  subroutine  input  data  and  output  are  described  and  samples 
are  given,  together  with  a  listing  of  the  subroutine. 

1.  Insertion  Procedure 

A  wave  propagation  code  normally  has  four  main  categories  of  operations: 
reading  the  input  data,  initializing  a  finite  difference  grid,  performing 
calculations  for  each  time  increment  at  each  grid  point,  and  printing 
the  computed  Information.  The  subroutine  SHEAR2  is  written  to  correspond 
with  this  organization;  it  has  parts  for  reading  data  and  initializing, 
for  calculating  stresses  and  damage,  and  for  printing  the  results.  The 
part  of  SHEAR2  to  be  used  at  each  CALL  is  indicated  by  the  parameter 
NCALL.  NCALL  can  have  the  following  values: 
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0  Read  data,  initialize,  and  choose  the  weighting  distribution 
of  crack  orientations  in  the  solid  angles  according  to  the 
value  of  NANG,  which  is  read  in  the  BFR  array, 

1  Read  data  and  initialize  only  for  one  material. 

2  Calculate  only. 

3  Calculate  and  print. 

4  Print  only. 

For  each  material,  SHEAR2  must  be  called  once  to  read  data  and 

initialize  variables.  For  the  first  of  these  calls,  NCALL  =  0;  for 

> 

subsequent  ones,  NCALL  =  1.  This  CALL  could  appear  anywhere  in  the 
initialization  phase  of  the  calculations;  it  occurs  in  INPUT3  of  TOODY 
and  in  GENRAT  of  PUFF.  These  CALLS  are  placed  so  that  the  SHEAR 2  data 
appear  with  the  other  data  for  the  material.  SHEAR2  can  then  be  called 
for  stress  calculations  alone,  printing  alone,  or  calculations  and  printing 
together.  These  latter  CALLs  should  appear  after  the  strain  increment 
has  been  determined  for  each  cell.  In  T00DY3,  the  second  call  is  for 
calculation  and  printing  and  is  within  the  loops  over  each  cell  and  each 
time  Increment  in  the  main  program.  In  PUFF  the  computing  and  printing 
CALL  is  in  HSTRESS,  the  subroutine  that  controls  all  stress  calculations. 
There  NCAIJ,  is  normally  set  to  2,  but  on  print  cycles,  NCALL  is  set  to 
3.  The  print  CALL  is  made  in  conjunction  with  the  computing  CALL  so 
that  certain  temporary  variables  that  are  computed  can  also  be  printed 
without  the  necessity  of  storing  them. 

The  formal  parameters  of  the  subroutine  are  listed  below  in  the 
order  that  they  appear  in  the  CALL  statement: 

SUBROUTINE  SHEAR2  (NCALL,  IN,  M,  K,  J,  IH3,  SX,  SY,  SZ,  SXY,  P,  TAU, 

DH,  DOLD,  DTO,  EH,  EOLD,  EN,  EMELT,  EP,  EX,  EY,  EZ,  EXY,  F,  YHL,  PLEN, 

ROT,  DROT,  ESC,  CN) . 


92 


» 

( 


Only  the  first  three  parameters  are  required  for  the  calls  in  which 
NCALL  ■  0  or  1.  The  meaning  of  the  parameters  is  given  in  Tables  3 
and  4.  Some  parameters  are  further  described  in  Section  4  of  this 
Appendix. 

2.  Outline  of  the  Subroutine 

The  subroutine  is  written  in  three  parts  for  reading,  computing 
stresses  and  damage,  and  printing.  In  the  reading  portion,  the  BFR 
and  NSIZE  arrays  are  read,  NSIZT  and  the  two  sign  factors  SS  and  SSE 
are  set,  and  the  weighting  factors  FNUC  are  initialized. 

In  the  second  part,  the  stresses  and  damage  are  computed.  First, 

the  T  array  is  constructed  from  the  existing  damage  quantities.  The 
z 

T  values  are  used  to  transform  the  external  stresses  and  strains  to 
z 

internal  values.  Then  new  stresses  are  computed  on  an  elastic  basis. 

If  yielding  occurs,  the  equivalent  plastic  strain,  plastic  work,  and 
plastic  shear  strains  in  each  band  orientation  are  computed.  Then 
existing  bands  are  allowed  to  grow  in  proportion  to  the  plastic  shear 
strain.  If  some  strain  remains  after  growth,  nucleation  of  new  bands 
is  computed.  Then  is  recomputed.  Finally  the  internal  stresses 
are  transformed  back  to  the  external  sign  convention  and  orientation. 

A  special  zero-stress  route  is  provided  if  fragmentation  has  occurred. 

The  final  part  of  the  subroutine  contains  the  printing.  Printing 
occurs  only  if  there  has  been  some  nucleation  of  shear  bands  in  the  cell. 
This  output  is  described  in  Section  6  of  the  Appendix. 

3.  Descriptions  of  Parameters 

All  the  formal  parameters  and  input  quantities  are  defined  in  Tables 
3  and  4.  The  subroutine  calculations  and  variable  names  closely  follow 
the  derivation  in  Chapter  IV  so  that  internal  variables  should  be 
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Table  3 


NCALL 

IN 

M 

K 

J 

IH3 

SX 

SY 

SZ 

SXY 

P 

TAU 

DH 

DO  ID 

DTO 

EH 

BOLD 

EN 

EMELT 

EP 

EX 

EY 

EZ 


DEFINITION  OF  FORMAL  PARAMETERS  IN  SHEAR 2 

An  Indicator  for  the  type  of  calculations  required; 
see  discussion  in  part  2  of  this  Appendix 

Data  input  file  number  (usually  5) 

Material  number 

Cell  number  in  the  K  direction 

Cell  number  in  the  J  direction 

A  damage  indicator  (20*TAU+2.9) ,  output  only;  initialized  at  2. 
Set  to  25  at  full  separation. 

2 

Deviator  stress  in  the  X  direction  (dyn/cm  ) 

2 

Deviator  stress  in  the  Y  direction  (dyn/cm  ) 

2 

Deviator  stress  in  the  Z  direction  (dyn/cm  ) 

2 

Shear  stress  (dyn/cm  ) 

2 

Pressure  calculated  for  the  cell  (dyn/cm  ) 

_  3 

Damage  (£  NR  )  in  the  cell  calculated  on  the  last  cycle; 
ranges  from  0  to  1 

3 

New  incoming  density  (g/cm  ) 

3 

Cell  density  determined  on  the  last  cycle  (g/cm  ) 

Time  step  (sec) 

Internal  energy  of  the  cell  in  the  present  cycle  (erg/g) 

Energy  of  cell  in  previous  cycle  (erg/g) 

Total  number  of  shear  bands  in  the  cell  (number/cm  ) 

Value  of  internal  energy  at  which  material  will  melt 
Total  effective  plastic  strain 
Strain  increment  in  the  X  direction 
Strain  Increment  in  the  Y  direction 
Strain  increment  in  the  Z  direction 
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EXY 

F 

YHL 

PLEN 

ROT 

DROT 

ESC 

ESC(M,  1) 

ESC(M,2) 

ESC  (M ,  3)  , 
ESC (M,4) 

ESC (M,5) 
ESC  (M ,  9) 
CN 


Tensor  shear  strain  increment 

Thermal  softening  factor;  varies  from  0  to  1 

Yield  value  calculated  on  last  cycle  (remains  unchanged 
(1)  if  the  yield  has  not  been  exceeded,  or  (2)  if  it  has 
been  exceeded  and  work  hardening  is  zero)  (dyn/cm  ) 

Work  devoted  to  plastic  strain  on  previous  cycle  (erg/g) 

Total  angle  of  cell  rotation  as  of  the  previous  cycle 

Additional  angle  of  cell  rotation  since  the  last  cycle 

Equation  of  state  variables  (array  (6,20)  must  be  dimen¬ 
sioned  in  the  calling  routine;  some  of  the  elements  of  this 
array  are  not  used  in  SHEAR2) 

Original  density  (g/cc) 

2 

Bulk  modulus,  C(dyn/cm  ) 


D  and  S  in  the  pressure  equation: 

P  ■  C/1  +  D +  S fX  ,  where  fi  =  DH/ESC(M,1)  -  1. 

2 

Shear  modulus  (dyn/cm  ) 

»• 

Gruneisen's  ratio. 

Array  contains  numbers  of  shear  bands,  band  lengths,  and 
other  damage  factors  for  each  cell.  The  number  of  array 
elements  required  for  each  cell  is  2p  NSIZE(i)  +  NANG; 
array  must  be  dimensioned  in  the  calling  routine. 
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Table  4 


BFR  array 
BFR (M,22) 
BFR  (M,  23) 
BFR (M, 24) 

BFR(M,2S) 

BFR (M, 26) 
BFR (M, 27) 
BFR (M, 28) 
BFR (M, 29) 
BFR (M, 30) 
BFR  (M,  31) 
BFR  (M,  32) 
BFR (M, 33) 
BFR (M, 34) 
BFR (M, 35) 

NSIZE  array 


INPUT  DATA  FOR  SHEAR2 

Nucleation  and  growth  parameters  for  shear  bands 
Growth  coefficient,  C 

G 

Growth  threshold,  e  (unused) 
go 

Nucleation  size  parameter,  R  in  the  relation  N  =  N 

exp(-R/R  )  (cm)  n  go 

n 

2  3 

Nucleation  rate  coefficient,  C  (sec  /cm  ) 

N 

Nucleation  threshold,  e 

no 

Ratio  of  displacement  to  band  radius,  b 

Maximum  nucleation  band  size  (cm) 

Ratio  of  successive  intervals  between  band  sizes 

Time  constant  for  deviator  stress  relation,  T  (sec) 

2 

Work  hardening  modulus,  Yd  (dyn/cm  ) 

Number  of  orientations  used,  NANG 

2 

Thermal  conductivity  (erg/sec/cm  )  (unused) 

Specific  heat  at  constant  volume  (erg/g/°K)  (unused) 

Critical  internal  energy  required  in  the  shear  banding 

material  to  permit  motion  of  the  band,  E  (erg/g) 

cr 

Number  of  size  groups  for  each  orientation 
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Identifiable.  However,  some  additional  description  is  given  here  of 

the  paraaeters  involved  in  the  size  groups  and  of  the  nucleation  threshold. 

BFR(M,32)  =  NANG  is  the  number  of  orientations  that  may  contain 
shear  bands.  We  have  chosen  six  orientations  for  the  present  cylindrically 
symmetric  and  axisymmetric  problems.  In  same  problems  fewer  could  be 
chosen  if  some  orientations  would  not  be  activated  because  of  symmetry. 

The  six  chosen  are  for  bands  whose  planes  are  normal  to  the  radial, 
circumferential,  and  axial  directions,  and  at  45°  between  radial  and 
circumferential,  between  radial  and  axial,  and  between  circumferential 
and  axial.  Since  only  the  latter  three  are  active  in  the  cylindrically 
symmetric  cases,  NSIZE  was  set  to  zero  for  the  first  three  orientations. 

The  direction  of  most  interest  is  at  45°  between  radial  and  circumferential. 
Probably  satisfactory  results  would  be  obtained  for  the  cylinders  if 
only  that  orientation  were  provided.  Currently  only  six  orientations 
are  fully  implemented  in  the  subroutine.  More  are  probably  not  needed 
until  complete  3D  problems  are  attempted. 

In  each  orientation  a  number  of  size  groups  is  provided.  The 

number  of  size  groups,  NSIZE,  should  be  large  (say  8  or  10)  if  considerable 

detail  is  desired  in  the  shear  band  size  distribution  for  that  orientation; 

however,  NSIZE  could  be  2  or  3  for  less  important  orientations.  For 

all  shear  bands  in  a  material  there  is  a  maximum  band  size  for  nucleation, 

BFR(M,28).  The  choice  of  maximum  band  size  should  be  larger  than  needed 

so  that  the  results  do  not  significantly  depend  on  this  choice.  In  a 

3 

shear  band  calculation  if  NR  from  the  largest  size  group  is  a  small 
fraction  of  ENR3,  then  the  computed  results  are  not  dependent  on  the 
largest  radius  size  selected. 

The  radii  of  the  bands  are  initialized  so  that  the  ratio  of  successive 

intervals  between  radii  is  equal  to  BFR(M,29)  ■  X,  an  input  quantity. 

Then  if  the  maximum  radius  is  R T  ■  BFR(M,28),  the  (N-i)th  radius  is 

N 
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R  -  R  (1  -  XN"S/(1  -  xV 

N—  1  N 


The  values  of  NSIZE  and  BFR(M,29)  should  be  coordinated  to  give  an  accurate 
treatment  of  the  band  size  distribution.  A  small  ratio  of  say  1.2  to  1.6 
could  be  used  with  large  values  of  NSIZE.  A  small  NSIZE  would  require 
a  large  BFR(M,29)  to  span  the  range  of  radii. 

The  threshold  plastic  shear  strain  is  BFR(M,26)f  the  value  of 
strain  that  must  occur  in  each  orientation  before  nucleation  can  begin 
there.  The  threshold  does  not  apply  to  the  equivalent  plastic  strain, 

—  D 

e  ,  but  it  applies  to  each  element  of  the  tensor  eP^  -  We  note  that 

( p  H 

ePs  is  only  the  shearing  portion  of  the  deviatoric  plastic  strain. 

<P% 

This  shearing  strain  is  obtained  from  the  flow  rule  or  proportionality 
between  plastic  strain  rate  and  deviator  stress. 


d,.j  ’  V* 

By  squaring  this  equation  and  summing  contributions  from  all  the  ij 
combinations,  we  obtain 


-  “*>2 

But  these  soma  are  proportional  to  the  squares  of  the  effective  stress 

—  -p 

and  strain  quantities  ((T  and  dC  )  defined  in  Equations  4  and  5. 

Hence 

dA  =  d€P/a 


so  the  flow  rule  can  be  written 

/ 


"ij  ■ df 


P  °1J 


or 


-ps 

d?iJ 


where  O  is  the  equivalent  or  von  Mises  shear  stress  and  (T. 

is  the  shearing  component  of  the  deviator  stress.  This  shear  component 
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is  found  by  taking  the  vector  difference  between  the  total  deviator  S 
and  the  normal  component  S  : 


2 

n 


In  general,  S  and  S  are  defined  as 


S  =  AS  +mS  +  n  S 
n  x  y  z 


S  =  I  O  +  nO  +nOr 
x  x  yx  xz 


S  =  I  ff  +  m  0  +nO 
y  xy  y  yz 


S  =  l  a  +  mo  +  nCT 
z  xz  yz  z 


where  O  ,  o  ,  O  ,  O  ,  O  ,  (X  are  the  stress  tensor  components  and 
x  y  z  xy  xz  yz 

I,  i,  n  are  direction  cosines  to  the  plane  of  interest.  S  ,  S  ,  S  are 

x  y  z 

x,  y,  z  components  of  the  stress  acting  on  the  plane.  A  standard 
mechanics  text  should  be  referred  to  for  details  of  this  transformation. 

As  a  guide  in  selecting  values  of  BFR(M,26)  based  on  measured 

thresholds  for  Sp,  consider  the  case  where  de  =  -  de,  and  de„  =  0. 

2  13 

This  strain  configuration  is  approximated  in  the  cylindrically  symmetric 

o 

calculations.  For  this  case  the  shear  strain  at  45  between  the  directions 

1  and  2  is  de?!!  =  de*3/  »/3  .  Hence,  for  an  observed  threshold  of  ep  =  0.28, 
12  v 

the  appropriate  value  of  BFR(M,26)  is  0.28/ ^3  =  0.16. 

The  major  new  material  parameters  are  C  ,  R  ,  C  ,  e  and  b;  each 

G  n  n  no 

represents  a  well-defined  physical  process  or  quantity.  From  the  present 
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computed  results,  It  appear?  that  C  Is  the  most  Important  for  determining 

n 

the  shape  of  the  fragment  size  distribution. 

4.  Transformations 

As  they  enter  the  subroutine,  all  stress  and  strain  quantities 
undergo  two  transformations — one  in  orientation  and  one  in  sign.  The 
orientation  transformation  allows  all  computations  in  the  subroutine  to 
occur  in  a  coordinate  system  that  rotates  with  the  material.  (The  shear 
band  size  groups  and  the  cumulative  plastic  strain  tensor  also  rotate 
with  the  material.)  The  sign  transformation  occurs  so  that  the  internal 
sign  convention  need  not  match  the  sign  convention  of  the  calling  program. 

Material  rotation  is  provided  by  the  calling  program  with  two 

variables:  ^  and  (positive  counterclockwise) .  The  ^  is  cumulative 

and  must  be  stored  external  to  the  subroutine.  Transformations  use  the 

angle  =  0  +  tiP/2.  Stresses  are  transformed  from  O  ,  (J  ,  O  ,  O  to 

x  y  z  xy 

the  starred  quantities  inside  the  subroutine  by  the  following  equations: 


* 

a 

X 


a  +  a 

X _ l 

2 


cos  2^/  +  O  sin 
xy 


(36) 


* 

a 

y 


a  +  a 

_X _ £ 

2 


cr  -  a 
_x _ y 

2 


cos 


2 V 


O  sin  2  ^  ' 

xy 


(37) 


a*  =  a  (38) 

Z  Z 


*  or  -  a 

a  =  -  — — - — 2L  sin  +  a  cos  2 
xy  2  xy 


At  the  end  of  the  calculation  the  internal  stresses  are  rotated  back  to 
the  external  configuration.  The  equations  for  this  final  transformation 
are  the  same  as  Eqs.  36  to  38  with  the  starred  and  unstarred  quantities 
interchanged  and  $'  changed  to  . 
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The  strain  increments  are  also  rotated  to  the  material  directions, 
and  the  same  transformation  equations  are  used  for  strain  as  for  stress. 

The  sign  transformation  is  made  to  account  for  a  different  sign 
convention.  The  internal  convention  is  positive  in  compression  for 
stress,  pressure,  and  strain.  If  the  convention  matches  the  external 
one,  SS  and  SSE  (for  stress  and  strain,  respectively)  are  set  to  +1. 

If  the  external  convention  is  that  tension  is  positive,  then  the 
appropriate  factor,  SS  or  SSE  is  set  to  -1  in  the  DATA  statement  in 
SHEAR2.  It  is  assumed  that  pressure  is  always  positive  in  compression. 


5.  Output  of  the  Subroutine 

A  sample  of  the  output  from  the  subroutine  for  four  cells  is  given 
in  Figure  40.  This  output  is  from  a  PUFF  run  with  cylindrical  symmetry 
in  which  orientations  4,  5,  and  6  are  activated.  The  parameters  have 
the  following  meanings. 


CN  Cumulative  number  of  bands  in  an  orientation. 

The  bands  at  each  radius  are  stored  and  used 
in  computations,  but  the  cumulative  values 
are  printed  for  ease  in  examining  fragment  size 
distributions. 


CL 

TAUZ 

TAU 

EP 


Shear  band  radius 

£  NR  in  each  orientation,  summed  over  size  groups 

£  TAUZ  summed  over  all  orientations 

£Aep  imposed  on  the  material,  summed  over  time. 


Cumulative 
plastic  strain 

EN 

K,  J 


£AeP|  in  each  orientation,  summed  over  time. 
Total  number  of  shear  bands  in  all  orientations. 
Coordinate  designators,  used  for  printout  only. 


101 


CM  « 

CM  ^ 

«  d 

CM~ 

CM  O 

n  o 

(M  m 

m  pH 

m  o 

m  ph 

m  o 

o  o 

O  O 

oq 

o  a 

1  o  o 

o  o 

o  o 

O  O 

o  p 

O  O 

o  o 

•  * 

•  ♦ 

•  «  *PH 

•  4 

1  ft  ♦ 

•  ♦  IO 

•  * 

•  « 

l  ♦  <o 

1  ft 

•  ♦ 

UJ  Ui 

UI  UI 

uiui  o 

UIUI  UI  ui 

UIUI  o 

UI  u} 

uj  ui 

UIUI  o 

UI  UI 

0  0- 

O  CM 

♦ 

m  o  mm 

KO  ♦ 

CM  CM 

m  © 

CM  CM  ft- 

CM  K 

K  O 

ui 

0  pi 

0  n 

OK  Ui 

cm  o 

cn  0  ui 

©  O 

0  in  ui 

©  pH 

O  P« 

PCM 

00 

i M  ru 

mo  in 

CM  CM  p-  K 

0m  m 

*-•  CM 

0  p-< 

CM  ©  0 

o  m 

m  in 

CM  CM 

•  • 

•  •  Q 

•  4 

1  •  • 

•  •  K 

•  • 

0 

• 

I  • 

0 

• 

CD 

I  CM 

O 

© 

o 

,  ® 

© 

m 

CM 

N  H 

NON  ■ 

pH  H 

i  CM  o 

CM  O  CM  « 

CM  H 

CM  «0 

CM  O  CM  ft 

CM  ph 

CM  O 

o  o 

O  O 

o  o  •  a 

o  a 

1  o  o 

o  ©  •  a 

O  Q 

o  o 

o  a  •  a 

O  O 

B  ♦ 

•  ♦ 

t  *i  UJ 

•  « 

1  ft  « 

•  ♦  UJ 

B  1 

•  ♦ 

•  ft  uj 

1  ft 

f  ft 

UI  UJ 

Ui  Ui 

UI  UI 

UJUJ 

|  UJUJ 

UIUI 

Ui  UJ 

ui  Ui 

UJ  UJ 

UI  Ui 

UJ  Ul 

pH 

in  k 

•© 0 

OK  p4 

K® 

l  on 

pH  0  IH 

•©  o 

in  cm 

OO  O 

o  in 

CM  0 

0  pH 

o 

o  ^ 

1  ®  o 

cm  ph  a 

0  « 

CM  o 

0  0  O 

O  o 

0  in 

m  in 

0  ^ 

CM  K  ♦ 

•H  (M 

i  in  © 

cm  m  ♦ 

K  CM 

CM  *4 

pH  0  ft 

0  CM 

cm  m 

cm 

•  ft 

•  «  Ui 

0  ®  K 

•  « 

•  • 

•  •  ui 

in  OD  CM 

•  ft 

ft  ft 

ft  ft  Vii 

•  • 

•  • 

o  m 

CM  ~  pH 

in  ph  o 

m  ph  rj 

m  o  ph 

m  o  -4 

m  o  k 

CM  O 

•  ©  * 

•  o  • 

•  o  ft 

•  O 

(MO  O 

ph  pH 

r-l  O 

ph  o  o 

M 

CVJ  cJ 

CM  O  O 

CM  O 

CM  O  O 

o  o 

o  o 

OO  C 

oq 

1  O  O 

0*0  o 

•a  o 

o  o 

o-o  o 

OO 

O  O 

o  o  e 

1  « 

1  ♦ 

•  ♦  •  N 

ft  « 

1  ft 

B  ft  •  M 

1  ft 

3 

UJ  UJ 

Uj  Ui 

UJUJ  3 

UIU| 

Ui  UI 

UIUI  3 

ui  uj 

Ul  Ui 

UI  Ui  3 

ui  Ui 

Ul  Uj 

IT  <d 

r<  0 

®5)  H 

m  0 

®  m 

«-  ®  <4 

•  K'O 

0  0 

i©  0  0 

p-  K 

C  K 

Op- 

(*)  h- 

o  ® 

®  i  *- 

©  0 

CD  K 

OO  K 

in  n 

'©  pH 

0  0  B— 

m  0 

©  m 

m  ° 

0  PH 

n  0 

0  ©  OK  ' 

m  ^ 

pH  0 

-•  CM  CM  K 

CM  CM 

O  0 

m  CM  0  K 

ph  CM 

0  CM 

•  ft 

•  i  o  m  m 

•  i 

•  ft 

•  4  CD  in  -h 

•  4 

•  •  o  «n  o 

•  • 

(m  o  a 

1  CM  O  O 

CM  O  O 

cm  m 

♦ 

•  o  « 

•  O  * 

1  o  UJ 

1 

O  UJ 

O  Ui 

o 

o  © 

. 

o 

o  0 

fl 

O  H 

rs 

O  •  CD 

PH  p4 

*h  o 

O  O  »H 

PH  Q 

pH  o  •  K 

•H  O 

pH  PH  ft 

e  d 

o  o 

O  O  pH 

O  o 

o  o 

e  o  ni 

o  o 

o  o 

O  O  0 

O  O 

o  o 

o  o 

>  • 

»  *  i 

•  ♦ 

•  ♦ 

♦  ♦  ft 

•  ♦ 

B  « 

B  ♦  • 

B  ft 

B  ft 

•  B 

UJ  ui 

UI  UI 

UIUI 

UI  u 

UIUI 

Ui  UJ 

UI  Ui 

UI  UJ 

UJ  UI 

UIUI 

UJUJ 

4  fl 

k  « 

CM  ©  ©  K  | 

K  0 

0  2 

OO  O  K 

0  K 

0  o 

O  0  O  K 

O  PH 

©  0 

00  o  in 

-•  m 

^  i*i 

<*>  »  o  &  n| 

-*  i| 

*  ai 

O  (VI  O  GD  ft, 

CM  O 

pH  C*\ 

—  ®  O  O  ■ 

CM  pH 

m  in 

o  in  o  0 

PH 

N  0 

in-0  o  0  2j 

0  p4 

0  m 

HlVJ  0^7 

©  ph 

CM  O 

CM  pft  O  0  Z 

m  cm 

CM  -H 

PH©  O  0 

.  .  .  ou| 

•  i 

•  ♦ 

•  •  «C  UI 

•  • 

■  ft 

•  ♦  re  UJ 

•  • 

•  • 

•  •  •  in 

o®  i 

O  0 

O  0 

o  m 

©  j 

© 

© 

K 

• 

1 

• 

ft 

• 

O  PH 

o  o 

o  p 

O  PH 

^  p 

1  PH  O 

O  pH 

PH  O 

PH  O 

p  pH 

O  O 

o  o 

o  a 

o  o 

o  a 

>  e  a 

O  O 

o  o 

O  o 

o  O 

O  © 

©  © 

♦  ♦ 

o  o 

ft  ft 

i  ♦ 

ft  ft'  o  o 

1  ft 

B  ♦  O  O 

•  ft 

B  B 

1  B 

©  c 

Ul  Ui 

UIUI 

UIUI  oo 

UJUJ 

UJ  UJ 

uiui  oo 

uj  ui 

Ul  Ui 

UIUI  O  o 

UIUI 

UIUI 

UIUI 

o  o 

CM  pH 

o  p 

OK 

o  o  • 

©cvj 

o  n 

ITO  oo  * 

m  0 

O  0  Oo  • 

K  CM 

CM  © 

CM  0 

©  o 

0  CM 

0  0 

OP 

•  o>® 

O  CM 

00  •  0*0 

CM  0 

0  m 

K’PH  ft  ©<© 

©  K 

O  0 

CM'© 

•  o 

PH  p 

cm  m 

oo 

PH  H 

0  n 

1  pH  PH  ©  © 

PH  PH 

0  K 

0  ph  O  O 

©  PH 

o  0 

CM  0 

O  o 

•  • 

o  * 

ft  ft 

•  4 

ft  •  O  M 

•  ft 

•  f 

•  •  o  H 

•  ft 

•  • 

•  • 

o 

O 

ft  K 

o  O 

ft  K 

o  o 

oS 

ft 

o 

ft 

a 

a 

IT 

o  O 

o  o 

PH  O 

oo  o 

O  p 

e  o 

cm  ©  oom 

O  PH 

O  O 

o  PH  o  OK 

O  PH 

pH 

—  — 

©  c 

o  o 

o  o 

OO 

OO  CM 

■O'O 

o-o 

0*0  oo  (VJ 

■  ©  e 

o  o 

•  O'O’O  O  PH 

0*0 

*0-0 

o-e 

©  C 

ft  ft 

o  o 

ft  ft 

ft  ft 

♦  ♦  o  o 

♦  ♦ 

ft  ft 

ft  BOO 

ft  ft 

B  1 

«  1 

©  © 

Ul  Ui 

UJ  Ui 

UJUJ 

ft  © 

UJUJ 

UJ  Ul 

UIUI  ft  o 

ui  ui 

Ul  Ui 

UIUI  •  o 

Ul  Ul 

UIUI 

oo 

K»0 

le-Kioo 

©CM 

0  © 

cm  moo 

K  IT 

®«PH 

o  ©•©•© 

•  ©  m 

m  ph 

f»lp 

ir  a> 

cm  m 

•e  1 

CM  CM 

*C  © 

ICO  '©  M 

-pie 

*C  CM 

0  O  <C  K 

m  cm 

CM  CM 

OCV 

0  o 

CM  0 

pH  CM 

•  m 

m  o 

0  CM 

cm  p®  •  m 

CM  ph 

pH  © 

m  ©  •  m 

0  in 

m  cm 

• 

ft  ft 

•  • 

•  • 

o  X 

ft  ft 

•  ft 

•  •  OX 

•  ft 

ft  ft 

•  •  ox 

•  • 

•  • 

•  • 

© 

FIGURE  40  SAMPLE  OUTPUT  FROM  SHEAR2  FOR  FOUR  CELLS  CONTAINING  SHEAR  BAND  DAMAGE 


6 .  Input  Data 

Four  sets  of  input  data  are  provided  in  Figure  41  as  samples  of 
the  use  of  the  subroutine  SHEAR2.  All  were  used  with  SRI  PUFF8  in  the 
one-dimensional  cylindrical-symmetry  mode. 

In  Shot  5084-3,  the  explosive  pressure  is  simulated  by  an  exponentially 

decaying  pressure  on  the  inner  boundary.  The  expression  used  is 

P  =  P6*  exp(T/T6) ,  where  P  is  applied  pressure  and  T  is  time.  The 

parameters  P6  (initial  pressure)  and  T6  (time  constant)  are  given  in 

the  NAMELIST  statement  at  the  end  of  the  deck.  For  Shot  5084-4,  the 

low  density  PETN  explosive  is  modeled  by  a  pressure- volume  relation 

4 

given  in  tabular  form.  The  two  experiments  of  Crowe  contain  COMP  B 
explosive,  which  is  treated  by  a  polytropic  gas  law  following  a  simultan¬ 
eous  detonation.  Input  units  are  in  the  dyn,  cm,  g,  sec  CGS  system. 

Each  data  set  contains  three  parts.  First  is  the  general  running 

and  printing  information  including  JCYCS  and  TS,  the  number  of  computational 

cycles,  and  the  problem  time  at  which  the  calculation  is  halted.  In  the 

second  part,  the  data  for  each  material  is  given.  The  EQST  card  contains 
"  _ 

the  Mie-Gruneisen  quantities  C,  D,  E  ,1,  H,  and  S  where  C,  D,  and  S 

sub 

are  the  Hugoniot  coefficients;  E  is  the  sublimation  energy;  and  F 

sub 

and  H  are  Grunei sen's  ratio  in  the  solid  and  vapor  states.  The  SH  cards 
contain  the  14  BFR  numbers  in  order.  Yield  strength  and  shear  modulus 
are  on  the  YIELD  card. 

The  layout  of  the  materials  and  cells  is  handled  in  the  final  set 
of  cards:  NLAYER  is  the  number  of  layers,  and  JMAT  gives  the  material 
number  in  each  layer.  After  the  NLAYER  card  is  a  card  for  each  layer 
describing  the  size  and  number  of  cells  in  that  layer.  In  these  samples 
the  cell  sizes  are  uniform  in  each  layer. 

7.  Listing  of  SHEAR 2 

The  subroutine  SHEAR2  is  listed  in  Figure  42. 
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I  DENT  FRAS  AND  FOB  ERLICH-3 
NTEOT  ■  6  N JED IT  ■  2 

JEDITS  *  1  3-5  7  9  U 

JEOIT  2  S2.P.COM3,  I  .1 1 1 .3*1 ,5» 1 .7 
NEOIT  «  in  JCVCS  «  300 

NMTRLS  ■  .4  MATFL  *  -2 

HF-1  SH  BAND  RHOS  ■  7.85E0 

EOST«  l.S89E»12  5.170E*12  7.360E*l0 

SH2  3.000E*01  »2000Etw0  l.IOOE-Oi 

1.4  3.000E-0H  3.000E*08 

NSIZE  0  6  0  H  A  8 

YIELD  *  1 .030E»l6  B . 190E* 1 1 

BELT  >  1 .o8SE*1 n  1.000E*08  0.E0 

PMMA-RKR  (BANKER)  RHOS  ■  1.184E*00 

EOST  ■  7.000E*10--4.050Etil  1.000E*t0 

YIELD  ■  1.000£«0ft  1.950E«10  2.850E*09 

HF-1  RHOS  ■  7.85E0 

EQST*  1.589L»12  5a70E*12  7.3fcQE*10 

YIELD  >  1.030EO0  8.190E*ll 

MELT  c  1.08SE*ln  1.000E*QB  O.EO 

LEAD.  (KOHN)  RHOS  ■  ...  11.355 

EOST  ■  5.o0nE*11  4.986E*1 1  9.155E*09 

NLAVERS  «  5  JMAT  *  «  1 

NZOnES*  1.  i  CELLS  Lb  1.512E*00 

NZONES*  1  5  CELLS  IN  8.750E-01 

NZONES>  1  5  CELLS  IN  0.508E«00 

NZONES*  1  in  CELLS  I*L.  &.500E400 

NZONFS*  1  5  CEILS  IN  1.000E*00 

EXTRA 

SNL1ST  P6  ( L)  *81  .5E9.T6  ( 1  ).■“!  .36E-5S 


FIGURE  41  INPUT  DATA  FOR  4  CYLINDRICAL 
INCLUDING  SHEAR  BANDING 


NALPH#  «  2 

12  la  14  15  IB  22  27  3o 

CKS  «  3.000F *01  TS  «  B.OOOE-i'S 

UZERD  «.  0.  _  _  NSCRb  *  OOOOll 1111 

CFP*  030  OPY*  002  NVAR  >  58  NCONs  0 

I. 69E0  0.25E0  5.17CE*13 

J. OOOE-03  C.17  0.C7O £*00  O.07r£*0O 

6.  ,?  .17  7.oonE*OV 

0  0  0 

9.975E-01 l.OOOF-Ol 
CFP  ■  000  DPV  «  001 
1.0DO£«(IO  2.500E-01  3964cE«ll 

CFP*  000  OPV*  002  NCGIv=  l> 

1.69E0  .  JU25E0  5.170E*13 

9.975E-01  l.OOOF-Ol 
CFP  •  000 -DP Y  *  000 
2.2  0.25  ?.01VE*12 

2  3  4 

INCH  ..  . 

INCH 

INCH 

INCH . - 

INCH 


PUFF  CALCULATIONS  OF  FRAGMENTING  ROUNDS 
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1DENT  FMI  WHO  FOR  ERLICH-A 

NTEDl  ■  p  NJfPIT  ■  l  N ALPHA  • 

JCD1I  •  SI «s2.COM3t  Mt».SiMi2.T 

NEOIT  •  1  fib  JCVCS  ■  BOO  CKI  •  3.00flF*01  TS  • 

NNTRLS  ■  A  MATFL  >  I  UZEAO  ■  0.  NSCRH  • 

PETN  JRHOS*  1.099E*00  C7P«  000  OPY»  000 

COST  •  1,0197*11 

IHAX  *  21 


J.OO'E-O* 

ooornim 


u;r>n.n 


ITTiWWJiffl 


nu-a  ii 


ITF'CTI  ■)!  f  Ul 


.4007*00  ,U4E*11  .3907*00  .7207*10  ,316F»00  .469t*n 

.2417*00  .3097*1 0  .200E*00  .2071*10  .1597*00  .1407*10 


IFiJXJ  U 

mimi 


tj.ij  n'ti 


.1411*04  .iloi-oi  .C4St«00  ,HK*0l 
.3147-01  .1301*60  .2817-01  ,040E*00  .200E-01 


. 102E*04 
«711E*0A 


luyi 


HF-l  SH  BAND  RHOS  ■  *>R5E0  CFP»  0J0  UPY*  002  NVAR  •  5P  NCOMa  u 

EOST>  l,5t9E*l?  5.170E«12  7,3407*10  l.OOEO  0.2BF0  3.170705 

SH2  3. V0OE*01  i200QE*u0  1.100E-02  1.000E-03  0.17  0.p?ut*00  u.07  F*oo 

1.4  3.000C-0H  3.0007*00  0.  .2  .17  /.00(E*09 

NS17E  OPOOPOOOO 

YIELD  m  1,03 Oittfi  0.100EO1 

MELT  •  l.l84E*ln  1.0007*9*  fl.EO  9.97SE-01  1.000F-01 

PMHA-PKB  (HArtKER)  RH04  ■  1.1047*00  CFP  •  000  liPY  •  001 

COST  ■  7.UOOEOO  4.0SOE«4l  1.000E*10  1.000E«00  2,5007-01  3«64(<E*11 
YIELD  ■  1.000E*0«  1.950E*10  ?.BSoE»09 

HF-l  RHOS  ■  7.0SE0  C7P»  000  DPY>  002  NCONa  o 

EOST>  1.5S9EU2  &.170E*12  7.34B£*lft  l.OOEO ...  .  0.2SE0  5.170E*13 

YIELD  •  l.i»3*E*l6  0.190E*11 

NfLT  ■  1.0P5E«ln  1.0007*00  0.E0  9.9T5E-01  1. 0007-01 

LEAD  (KOHN1  RHOS  ■  11.355  £7P  -000DPY  >  POO 

EAST  m  5,0047*11  4.4#*C«1!  9.1S9C«09  7.2  0,?S  2.01*E*1? 

PLAYERS  •  5  JMAT  ■  12345 

n/ONEs>  1  3  CELLS  IN  l.SJ2E*00  INCH 

N?ON7s>  1  A  CELLS  IN  P.750E-0!  INCH 

N20N7S-  1  3  CELLS  IN  1,12SE*Q0  INCH 

NZONES*  1  3  CELLS  IN  2.5.40E*00  INCH 

NZONES>  1  7  CELLS  IN  1.000F*00  INCH 

EXTRA 

0NL1S1  U(2)«l. o»0. .1. 0t0.«1.0*0.i\. 0*0. «1. 00 


FIGURE  41  INPUT  DATA  FOR  4  CYLINDRICAL  PUFF  CALCULATIONS  OF  FRAGMENTING  ROUNDS 
INCLUDING  SHEAR  BANDING  (Continue) 


I0E*T  FR 

4  /H4AMTG  PND  OF  HF-1 

TO  simulate  crones 

TEST*  1  AND  ?. 

ntedt  ■ 

NJEP1T  a 

1 

NREZON  a 

0 

NALPHA  a 

2 

JED IT  a 

S1»S?,C0M3.».  _2*1«2«3, 

2. 5, 2,7,2 

.9 

NEDIT  a 

ion 

JCVCS  » 

300 

CFS  ■ 

30.  EO 

IS  a 

rt.no. 'E-os 

NHTRLS  ■ 

? 

MATFL  ■ 

1 

U7ER0  a 

O.EO 

NSCRh  a 

COHPB 

RHOS_» 

1.T2E0 

CFP.a_000 

OPT  a  012 

EQST  a 

l.to 

O.EO 

1.E0 

1.B41E0 

1.P41EP 

OEXPL  a 

4.469E*lo 

TENS  a 

-1.000E*09 

O.EQ  _ « 

1.E0 

MELT  a 

-l.to 

HF-1  SH 

BANO 

RHOS  ■ 

7.85E0 

CFPa  030 

DPT*  002 

NVAP  a  58 

NCONc  J 

EASTo 

1.S89E  *12 

5.170E.12 

7.360E*10 

1.69E0 

0.25E0 

5.1 7t  E*13 

SM2 

3.COOE*Ol 

.2O0CE«b0 

1.100E 

-02 

1.000E-03 

0.17 

0.07t>E*00 

0.07'IE*00 

1.4 

3.000E-0B 

3.000E*08 

6o 

.2 

.17 

7.00(»E*09 

NSIZE 

Jl  -_5 

_.0  .8 

_  8 

8 

0  0 

0 

YIELD  a 

1.030F.10 

8.190E.U 

MELT  a 

1.0B5EO0 

],000E*UP 

O.EO 

9.975E-01 

1  .OOOF-01 

NLAYERS- 

■  ..  2 

JMFT 

__  1 

2 

nzones* 

1  10 

CELLS  IN 

1.5 

INCH 

NZONE$a 

1  10 

CELLS  IN 

8.750E 

-01 

INCH 

I DENT  FR 

S  FHAQ  POUND  OF  ARMCO  IRON  TO  SIMULATE  CROWES  TESTS  3  AND  4 

NTEoT  a 

P  NJEOIT  a 

1 

NREZON  a 

0 

NALPHA  a 

2 

JEDITS  ■ 

S 1 «  S2  *S3  .COM3 —  2.  It 

2.2j  2.3, 

2.4,  2.5, 

2.6,  2.7* 

2.8,  2.9 

NEOIT  a 

100  JCYCS  * 

300 

CKS  a 

3.000E.01 

TS  a 

rt.CO  iF-05 

NMTRLS  a 

?  MATFL  a 

1 

UZERO  a 

0. 

COMPS 

_ RHOS .« 

1.72 

CFP  M  £00 

OPY.a  01^ 

EOST  a 

1.  0. 

1. 

1.841 

1.841 

0. 

0. 

OEXPL  a 

4.489E.10 

TENS  a 

-l.000E.O9  fl.  .  - 

1. 

_ 

MELT  a 

-1. 

ARMCO  SH 

BANK  MHOS  « 

7.85E0 

CFPa  030 

DPY«  00? 

NVAR  a  58 

NCONc  <i 

EQSTa 

1.589E.12  5»170£a12  Z.3«0Eflil 

1.09E0 

0.25E0 

5.17ut*13 

SH2 

3.000E.01  .2000E»00 

1.100E-02 

3.000E-04 

0.17 

9.O70E.O0 

0.07:'E*0U 

1.4  3.000E-08 

3.000E«08 

6* 

.2 

.17 

7 .  OO'lF  *09 

NSIZE 

0  ft.  0  ,  ,8. 

..  8 .  8 

0  .0 

0 

MELT  a 

1.085E.10  6.460E»09 

5.700E-01 

7.000E-01 

1.R50F-01 

TIELO  a 

2.000E.A9  8, 190E.il 

NLAYERS 

a  ?_JMAI  ■ - 

1  2 

NZONES* 

1  lit  CELLS  IN 

1.5 

INCH 

NZONES* 

1  10  CELLS  IN 

0.78 

INCH 

FIGURE  41  INPUT  DATA  FOR  4  CYLINDRICAL  PUFF  CALCULATIONS  OF  FRAGMENTING  ROUNDS 
INCLUDING  SHEAR  BANDING  (Conducted) 
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SUBROUTINE  »HEAR2INCALL*IN*N»ll»J*tH3*»X*SV*S/*»XY*F»TAU»0H*«0L0* 
lQlU,EP.lQL0.Eh.lNELT*|F.EX*IJ»El*MT|F*YML»PLEN*RQ?»0ROT*ESVr*CN! 

ROUTINE  FOR  CONFUfATION  Of  ITREIMB  <»1TH  RATE-OEPENDCN! 
YIELD  MODEL  FOR  OCVJaTORS  A  NO  NIC«MUNCUCN  fON  WtlllWt 
if  THtiESHOLO  FLABTIC  STRAIN  IS  REACHED.  SHEAR  BANOS  ARE 
NOCLEaTCU  A  NO  BRORN  IN  6  ORIENTATIONS. 

SV»  ST*.  SI  AR«  pf VI ATOMS  IN  (XTKRNAt,  »IBN  CONVENTION, 
f  IS  POSITIVE  IN  COMfRCSSlON,  INTERNAL  SION  CONVENTION  IS 
POSITIVE  IN  COMfRCSSlON  fOR  ALE  STRESS  ANO  STRAIN  QUANTITIES 
ST  IS  TOTAL  STRESS  Al.fRfVIOOS  HMf.  SI  IS  NCR  OEVXATOR. 

EX*  EV*  CZ.  EXT  ARC  STRAIN  INCREMENTS  IN  EXTERNAL  SION  CON, 

SS*  SEE  CMANOC  EXTERNAL  SIBN  CONVENTION  TO  INTERNAL  fON 
STRESS.  ANO  STRAINi  JESRCCTIVCLT,  . 

DIMENSION  BFR(OOB)  *NS1ZC(30«B)  *FNUCIB>  .TAUHn)  tEFR (3)  »VFM(b>  • 
..  IST.I4)  *£*(♦>. Bf<*>»T|FIS>,CSCC6,20)*NI»T  IB) 

I  *CN(IVS) *OCF(Ai *CLA(ISSt  *CNA(I00l tVNAA(B) 

EQUIVALENCE  (CNA.CLA) 

c 

DATA  Sf *SSE/*1 •**!•> 

NCttNCICL*! 

_  SO  TO  IlJ»,l0.1#«*10«J*00)NCl 
"10  MCAOI IN, 1002)  A1*A2* (BFR(M.l) *1>22*3S) 

PRINT  lu02*  A| «AZ, (DFR(M.I) • I >22 *30) 

1002  FORNAMcAS*TC10«I/10X*TC1B,2)  . 

READ C IN,  1003)  A1  «A2*  INSUE  (M* II  t  1*1 *0) 

PRINT  lOOT* A | ,A2*(NSlZC(MfI)*Ial.9) 

1003  FORMAT (2AS*1*IS) 

VNAX(H)sO, 

NSI2T(M|aNSIZEfM*n 

_DO  I*  1*2 *S ...  . 

14  NS12T<M)aNSlZT(M)*NSlZE(N,II 
VFH(M)a|. 

IF  (NCALL  .CO.  11  00  TO  OS 

NANRaOFH (Mt 32) 

DO  16  !■! ,3 

FNucM>».mm 

16  FNUC 1 1 *3) a.ZZZ222 

IF (NANO  -  A)  20*00*30 
20  fNUCl”. 333333 

IF  (NANb  .BE. 4)  FNUC I  *.25 
00  25  Ia| ,4 
ZS  FNUC(l)aFNUCI 

IF  (NANb  .£0,2)  FNUC (21  M.6B66667 
IF  (NANb  ,Ne.  S)  SO  TU  30 
FNUC (4) ] Z5 
FNUC(T)a.l2S 
FNUC (5) B0, 

. Ft)UC1<M«0. 

00  TO  40 

30  00  39  laT.NANB 

FNUC ( 1 ) ■. 1 1 11 1 1 
3B  FNUC (1*3) a. Ill 111 

40  CONTINUE 
65  RETURN 
C**M 

c 
c 
c 

c***< 


COMPUTE  STRESS  AND  UANASC 


SHEAMZ 

Z 

SMEAMZ 

3 

SHEAMZ 

4 

SHEANZ 

s 

SHEAMZ 

6 

SHtANZ 

7 

SHEAMZ 

b 

SHEAMZ 

9 

SHEAMZ 

10 

SHFAMZ 

11 

SHEAMZ 

1Z 

SHtANZ 

13 

SHEAMZ 

14 

SHEAN2 

IS 

SHEAMZ 

16 

SHEAnZ 

17 

SHEAMZ 

16 

SHtANZ 

19 

SHEANZ 

20 

SHEANZ 

21 

SHEANZ 

Zz 

SHEAMZ 

23 

SHEAMZ 

24 

SHEANZ 

ZS 

SHEANZ 

Zb 

SHE  AMZ 

27 

SHEAMZ 

26 

SHEANZ 

29 

SHEANZ 

30 

SHFAMZ 

31 

SHEAMZ 

32 

SHtANZ 

33 

SHEAMZ 

3* 

SHFAMZ 

3S 

SHEANZ 

3b 

SHEANZ 

37 

SHEANZ 

3b 

SHEANZ 

39 

SHEANZ 

40 

SHEAMZ 

41 

SHFAMZ 

42 

SHt  AMZ 

43 

SHEANZ 

44 

SHEANZ 

45 

SHEAMZ 

4b 

SHEAMZ 

47 

SHE  AMZ 

4b 

SHEANZ 

49 

SHEANZ 

50 

SHF  AMZ 

51 

SHFAMZ 

5Z 

SHEAMZ 

S3 

SHEAMZ 

54 

SHEANZ 

55 

SHfARZ 

5b 

SHFAMZ 

57 

SHEANZ 

5b 

SHEAMZ 

59 

SHEAMZ 

60 

SHE  AMZ 

frl 

SHEARZ 

bZ 

FIGURE  42  LISTING  OF  SUBROUTINE  SHEAR2 
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mm 


100  IF  I  IN)  ,«[,  25)  «0  TO  000 

if  <vmaaimi  ,e«,  p.)  vMA»(M)«sooT(esc(MtSi/esc(M«in 
c  COMPUTE  StflCSS  REDUCTION  FACTORS  TAUZ(I) 

TAUaO* 

JN*0 

...  _  DQ  110  NO *1.  NANO 
TAUZ(N«)ap, 

IF (NSlZt (M.NG)  *EQ*  0  .OR.  CNI JN*1 )  .EQ.  0.160  TO  110 
NS1ZEN*NSIZE<M.N6) 

DO  120  IbI.NSIZEM 
JNN«JN*2»i-l 

120  T  AUZ  I NO.)  nTAUZ  I NG ) .CN ( JNN )  *CN  I JNN*  1 )  **3 
TAUaTAU*TAU2(Nfi) 

110  JN»JN*2«Nsl2f (N.NO 

a 

INITIAL  TRANSFORMATION 


ADJUST  SIGNS. ROTATE  STRESS. TRANSFORM  TO  STRESS  IN  SULIOCST) 
RT«ROT*ROT*OROT 
EMUaOOLU/ESC (M. 1 ) -1 . 

PH*EMU* (ESC (M«?) *EMU* (tSC (M.3) .£MU*ESC IM.4) ) ) 

PS«RH* ( 1 .-ESC (M.9) *EMU/2, ) «OOLD«CSC ( M,9>  *E0LD 
If  IRS  .GT.  0.)  Rap-PsaTAU 
.  SAaCSX.SV)*SS/2.«P 
SORaSlN(RT) 

CORaCOS(RT) 

_  SBa((SA-SY)/2.*C0R«SXY*S0R>*SS 
00  140  Ial.4 
140  STIll-O. 

G2a?.#ESC  CM. SI 

ST  III ■  ISA«SR)/AMAXl (0.02. ( 1 I3.*TAUZ I1)»1.5*(TAU2(4) *TAU7  1 5) ) ) 
1  aVFR(H) ) ) 

ST  1 21  ■  (SA-SB1/AMAX1  (0.02.  ( 1 (3.*TAUZ (2)  *1 .5*1  T AUZ  (4) *TAU?  t«»l  >  I 
1  aVFR(H))) 

STI3»a(P-(SX.SY)*SSI/AMA*l(0.02.<l.-(3.*TAU2(3).l,5*(TAUZ(5>. 

1  TAUZIAM>*VFR(M))) 

ST  141 ■ ( ISY"SX) /2.*S0R»SXY*C0R>  *SS/ 1 1 •“ 1 1 ,5* IT AUZ 1 1 ) »TAU2 1?) 1*3, 

1  *TAU2 (41 | *VFR(M) ) 

RaJSI (1) »ST (21 aST (3) 1 /3» 

c  rotate  strains  to  band  orientations 

EAa(EX.EY)/?. 

EBa (EX-4 Y ) /2.acOR.EXYaSQR 
EBARaO, 4007a (OH-DOLO) / (DH.POLD) 

NSTERaSwRT ( ( ABS (EA ) ♦ ABS CEO  I  * AOS ( EbAR) ) / . 002) 

NSUR«MAXo  INSTEP*  1  > 

ESIlia  (Ea»EB)*SSE/NSTEP 
ESI2I •  (Ea-EB) aSSC/NSTlR 
_  ...  E4 L31 43. #EBAR/NSTER“ES ( 1  leES  12 1  - 

ESI4) a  I ( IEY-EX ) /2.*S0R»EX Y»COR) aSSE) /NSTEP 
00  000  NSal.NSTER 
DO  100  lal .3 

100  SE I I ) a$T I I ) *02* (ESIZI -EBAR/NSTER ) -R 
SC  14) •  ST (4) »Q2*ES<4) 

SNaSQRT (l ,5*(SC(1 )*#2.SE (2J  **2.SE (3) *a?.2,*SE (4) **2) ) 

OHNaDOLO .FLOAT (NS) /FLOAT (NSTEP) * (DH-DOLD) 

EMUaOMN/ESCIM.D-l. 

PHaEMUa (EpC (M.2) .EMU* (ESC (M.  3) .EMU*ESC (M.  4))) 

PE»RM*d.«ESC  (M.  0)  *EMU  /2.  )  *PMNaCSC  IM»  9)*EN 

YlayW. 

IF  (SN  ,LT.  VI)  60  TO  SOO 


SHEAR? 

t>3 

SHEAR2 

44 

SMEAR? 

f-S 

SHEAR? 

44 

SHEAR? 

47 

SHEAR? 

6H 

SHE  Ah? 

6*. 

SHF AH? 

70 

SHEAR? 

Tl 

SHEAR? 

7? 

SHF  AH? 

73 

Shf AH? 

74 

SHEAR? 

75 

SHEAR? 

74 

SHF  AR? 

77 

SHEAR? 

74 

SHEAR. 

79 

SHEAH? 

HO 

SHF  AH? 

HI 

SHEAR? 

H? 

SHEAR? 

H3 

SHEAR? 

84 

SHEAH? 

85 

SHEAR? 

6b 

SHE Ah? 

F  7 

SHEAR? 

88 

SHEAH? 

H9 

SHEAR? 

90 

SHEAR? 

91 

SHEAR? 

9? 

SHF  Ah? 

93 

SHf AR? 

94 

SHF  A 

95 

SHEAR? 

9b 

SHEAR? 

97 

SHEAH? 

98 

SHF Ah? 

99 

SHEAH? 

100 

SHE  AH? 

104 

SHEAR? 

It'S 

SHEAR? 

1  04 

SHEAR? 

107 

Shf  ah? 

108 

SHEAR? 

109 

SHFAR2 

110 

shear? 

111 

SHEAR? 

U? 

shear? 

113 

5HEAH2 

114 

SHEAR? 

US 

shear? 

116 

SHEANZ 

117 

SHEAR? 

1  18 

SHEAR? 

119 

SHEAR? 

1?0 

SHEAR? 

1?1 

SMEAR? 

1?? 

SHEAR? 

1?3 

SHEAR? 

l?4 

SHEAR? 

1  ?6 
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YIELD  *M0  PLASTIC  *T«*IN  CALCULATIONS 


FXPT  •  EXP ( -OTO/BFR |R»  30 1 /NBTEP) 

YEBa(Vl«BFR(M.31>*SN/?./02>/U.«WP<N.3t)/2./G2)/SN 

00  ISO  1.1,3 

ISO  SE(IlatXTl|>.P)*EXPT*<VCtMCm*S»M<N«SomMUI-STfl)*P>/ 
10T0*NSTbP|*(l.«CXPTI 

SC  (At  Ml (At  af XPT. ( yES*XC (At  *>FR (N.30I • (SC (At -ST (A) | /OTO*NsTCP) 
1  II. -CART) 

00  200  1.1,3 

200  OCP ( 1 1  afcS ( I > -CBAR/NSTEH- ( SC ( ! t -ST ( 1 1  *P t /#2 

OCP ( At aLS(A| • (SC (A) -ST (Alt /S2 

06AMMA.S0PT (1>S*(0CP( 1 1  **2*0PP(2>  *ap*0IP(3t*a2t  *0,7MDCP (At 
»••») 

YNLaVHL«8PP (N.31 ) aoOAMMA 

OPLCNRa ( (SC (It  *ST(1I-Pt  *DCPI I) • (SE (?) *ST (2)*P)*0CPt2l* (SE (3! *s 
lT(3t-P|aoeP(31*ISE(At«bT(A||aoEP(Alt/2*/OHNaAMARl(o.«l.-lAUt 

OPLCNNaatta  (DPLCNRI 

ep«cp*osamma 

PLCNaPLEN.OPLCNR 


COMMUTE  PLASTIC  strain  in  each  OSICNTATION 


STRlaABSISC (At  t 
STR2*AHS(SE(AI  I 
STR3ao. 

STRA  a AOS ( St ( 1 1 -SE ( 2 1  >  /? . 

STRSaSONT ( (SE ( 1 1 -SE (31 1 **2*2. »$E (a>  aa2t /2, 

STROaSQrtT ( (SE (2) -SE (31 l*#2*2.aSE (♦)••?) /?. 

SNasOHT ( 1 . 5* ( SF ( I I *«2 * SE ( 2 I **2 *SE ( 3 t **2 • 2 . *SE ( A I *#2 ) ) 

TCP ( 1 1 aUOANNA/SN«STRl 

TEP(2t  aLG»HHA/SN*STR2 

TEP(3>  aU0AMHA/SNa$TR3 

TEPtAlaUOAMHA/SNasTRA 

TEP(Sl  atiGAMHA/SN*STRS 

TEP(Alau0AMMA/SNaSTH6 


GROWTH  PROCESS 


NT0Ta2*NSl2T (Ht 
DO  2S0  lal.NTOT 
250  CNA ( 1 1 aCN ( 1 1 

IE  (£N  ,C<J.  01  GO  TO  3*0 

jNaO 

DCaVHAX (M) aoTO/NSYEP 
00  350  Nflal ,NANG 
DGAHao, 

IP  (NSUEIM.NGt  ,fO,  o.  .OR.  CN(JN*1)  .ED.  0  .OR.  T(P(NG>  .1 
I  GO  TO  JA5 

£XPCa£*P(BFR(M.22taTEP(N6>  t 
NSI2CNaNSIZC (N»NGt 
DO  300  |ai.*$!2EN 
jN2aJR*j(*  (NSIZE  (N*NG)*1"I) 

CL A ( JN2I aAHINI f CN ( JN2 1 *CXPE *CN ( JN2 1 t DC) 

300  OGAHbOGAM.cn ( JN2- t ) *3. 1**BER (M,27) • (CLA ( JN?) a*3-CN ( JN?) a*3) 
JP  (OGAM  ,lE.  TFPfNGDGO  TO  3AS 
RfcaTCMNGt/DOAM 
DCRaDCaHR 

rxPCaEXPfaaRR 

NSlZEMaNStZE (H.NG) 

00  3A0  la),NSIZEM 
JN2aJN-r*  (NSJ7F  (H*EG)*1“1) 


SHE  AH? 

I?*> 

SHI  AM? 

1*7 

SHE  AH? 

1?G 

SHE  AH? 

)?« 

SHE  A  HZ 

13o 

SHfe AH? 

131 

SH*  AH? 

1  i? 

SHE  AH? 

133 

•  SHEAH2 

13* 

SHFAM2 

1 35 

SHF AH? 

13b 

SHF AH? 

137 

SHF AH? 

138 

SHF AH? 

1  JO 

SHEAR? 

1*0 

SHE  AH2 

1*1 

SWAN? 

i»? 

SHE  AH? 

1*3 

SHEAR? 

1** 

SHEAR? 

1*5 

SHE AH? 

l*b 

SHF  AH? 

1*7 

SHE  AH? 

1*8 

SHEAR? 

1*0 

SHEAR? 

150 

she  AM? 

15) 

SHEAR? 

15? 

SHF AM? 

153 

SHF  ah? 

15* 

SHE  AH? 

155 

SHE  AH? 

15b 

SHE  AH? 

157 

SHF AH? 

158 

SHE AH? 

ISO 

SHEAR? 

160 

SHF AM? 

161 

shear? 

lb? 

SHEAR? 

183 

SHEAR? 

lb* 

shear? 

lb5 

SHE  AH? 

16b 

SHE  AH? 

167 

she  AH? 

166 

SHEA,*? 

189 

SHEAR? 

17C 

SHEAR? 

171 

SHEAR? 

i  n 

she  AH? 

1  73 

0.)  SHI  AH? 

17* 

SHEAR? 

175 

SHEAR? 

178 

SHE  AR? 

177 

SHEAR* 

176 

SHEAR? 

179 

SHE  AH? 

160 

SHE  AH? 

161 

She  ah? 

1**? 

shear? 

1*3 

she  ah? 

IP* 

SHF Ah? 

lb5 

SHE  AH? 

1H6 

SHEAR? 

167 

SHE  Ah? 

16* 
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340  CL  AC  JNZLpAMINI  (CM  I JN2.1  •£*  PE  tCN  I JN2I  OCR) 
3*5  TEP<N6l»ANAXl  (O.tTEPCNG)OGAMI 
390  jM«JM*N$I7t(M,M6)*2 

..  360  CONTINUE  . 

00  369  NGal (NANO 

365  CN1NT0T*NG) ■CM(MTOT  *NG)  *TEP(NG) 


MUCLCATIOM  PROCESS 
UPM*a 

DO  370  NGal .NANG 

170  TEPMaAMAXl (TEPMtCN (NT0T*NG) *TEP  t M6) I 
IF  iTIPH  .LI.  8FR<M(?6)I  00  TO  500 

JNM 

00  6S0  Neal .NANO 

IF  JMSUElM.NO)  .EO.  0)  00  TO  «5o 

IF  (CN(NT0T*N0)*TCP(N6)  ,LT.  BFR (M.26I  .OR.  TEPCNO)  .IT.  l.t-5) 

1  DN0*^CP  <n!>  OOP*  IM»|B)  PFNUC  CNG)  •tDPLINR/D i 06NSTEP2BFR  (M.3SI )  *•* 
CNPaO. 

NSIZEMaNSIZE IM.NG) 

00  440  lal.NSIZEM 
I IaNSIZt (MtNOI *1*1 
JNIaJN«?MI 

IFICLAUNI)  .ME.  O.i  GO  TO  420 

CLA I JNI I «BFR (M*28) * ( 1 .OFR  (M.  29)  I ) ✓ C 1 ,»BFR(M*29) #,NS1ZE (M.NG) 1 

CNC JNI). CL AC JNI) 

k20  CNL»ONO*EXP (• CCLA ( JNI 1 *CN ( JNI ) I /2./BFR (Hi24) ) 

JNN»JN.2*||-1 

CNA(JNN)aCNL-CNR*CN(JNN) 

*40  CNRaCNL 
ENa£N*CNL 

*50  jNaNSlZE (M«NO) *2* JN 
A70  CONTINUE 

CONFUTE  TAU  AND  REFILL  N*TN  ARRAYS 


TAUaO. 

JN-O 

IFIEN  .EO.  0.)  GO  TO  500 
DO  *90  NGal, NANG 
TAUZCNGtao. 

IF  (NSUElM.NG)  .EQ.O  .OR.  CNA(JN.l)  .EO.  0.)  GO  TO  *90 
NSlZEMaNStZE(MiNG) 

00  4«0  lal.NSIZEM 

JNNaJN.2*I-l 

CN ( JNN* U «CL A C JNN* 1 ) 

CN ( JNN ) ■ CN A ( JNN ) 

*GU  TAU2(NGI aTAUZ (NG1*CNA  C JNN1*CLA  < JNN*I)  ##3 
TAUaTAU«TAUZ(NG) 

*90  JNajN.N5IZE(M,NG)*2 

IF CTAU*VFR CM)  .GE.  1.)  GO  TO  000 
500  CONTINUE 
PaPE 

00  550  Ial.3 
550  ST (II «SE <  I  >  *P 
ST (4) *SE («l 
600  CONTINUE 

TRANSFORMATION  TO  GLOBAL  ORIENTATION 


SMEAR2 

1  HR 

SHE  AR2 

190 

SMEAR? 

19) 

SHtAR? 

192 

SHEAR? 

193 

SHE  AH 2 

19* 

SHE AH 2 

195 

SHE  *RZ 

196 

SHE  AH2 

197 

SHE AH2 

1 98 

SHE.AH2 

199 

SHE AH2 

200 

SHt AH2 

?01 

SHE  AH2 

?0? 

SHEAH2 

?03 

SHt AR2 

?0* 

SHEAR? 

2  OS 

SH»  AH2 

206 

SHEAR? 

?u7 

SHEAR? 

206 

SHEAR? 

?09 

SHEAR2 

210 

SHE  Ar2 

?ll 

Shear? 

?1? 

SHtAR? 

213 

SHEAR? 

21* 

SHEAH2 

215 

SHt AH2 

?  1 6 

SHE  Ah 2 

✓  17 

SHF  Ak2 

?1E 

SHEAR2 

?  1 9 

SHF  AH? 

?2t 

SHF.  AH? 

??1 

SHEAR? 

??2 

SHE  AR? 

??3 

SM>  AW? 

??4 

SHEAH2 

?2S 

SHEAR2 

??b 

shear? 

227 

SHfc  AH2 

??P 

SHEAR? 

??9 

SHF  AR? 

?3t 

SHE.AH2 

r3l 

SHtAR? 

?3< 

SHE Ar2 

23? 

SMEAR? 

23‘ 

SHEAR? 

?3‘ 

SHE  AR? 

?3t 

SHEAR? 

237 

SHF  An2 

?3K 

SHE  AR? 

?3'i 

SHEAR? 

?*( 

SHEAR? 

?*1 

Shear? 

?*? 

SMEAR? 

?*3 

SHE  »R2 

?*» 

SHE  AR? 

?45 

SHF AH 2 

?4b 

she:ar? 

2*7 

SHE An2 

?4il 

SHE.  AH2 

f  *7 
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STU)>S1U)MMAX!<0.tU.*ll.l*lT4UZm«UU2tt>>O.»T*U2t*>)»9F»(M>  SHEAR2 

?S? 

1  II 

SHFAH2 

253 

_ 

SI  111  "SI  dUANAXl 10  •  •  9*  *T4UZ  ill  *1«S*(TAU2  (41 tTAUZ  (511  l#VFF  (Ml  $HfAR2 

25* 

1  II 

SHEAR? 

?SS 

ST  (?) ■$! (2UANAX1 (0.* ( l«“ (3.*TAUZ (El *l»5*(TAUZ 141 *TAU? (Oil l#VFRIRI  SHE AR? 

?5r 

1  11  ...  .  . 

SHE  AM? 

?57 

ST (9I*$T (3I*AMAX1 10* » (!«• (3.*TAUZ (3) *1 .5*(TAUZ (SI •TAUZ(O) 1 )*¥FH(Hl  SHEAR? 

?5« 

1  II 

SNEAR2 

259 

•■1ST ( 1) *ST (t) *5J (311/2* 

SHEAR? 

260 

S*«(ST  < ll.ST (2) 1/2# 

SHEAR? 

261 

S8*<STU)-ST<2*)/2.»C0R-STt*)M0R 

SHEAR? 

?b? 

-  - _ 

SXVal  (SI  (li»$T(21 l/2.4S0N*ST  (019CORI  *8S 

SHEAR? 

?A3 

SX«(SX*Sl»Fl •tf 

shear? 

26* 

SV«(SA-kB«PlMS 

SHFAR? 

265 

_ _ 

S2«CST(9)-Pl«St 

SHEAR? 

?R6 

IF  (PE  .OT.  A.)  PbP*P£<*TAU 

SHEAR2 

?67 

!H3b?0.*TAU*?.9 

SHEAR? 

?6h 

IP  INCALL  .•€.  3)  00  TO  900 

SHEAR? 

270 

RETURN 

SHEAR? 

?71 

C#tM 

SHEAR? 

27? 

c 

CONPLETt  REPARATION 

SHFAR? 

?73 

£•••# 

•••••• 

SHEAR? 

?7* 

_*00 

ERUaON/'tsciH.n-i. 

SHEAR? 

?  75 

PHaCNU*  (ESC IH«Z) ♦EMU* (ESC (Mi  9) ♦EMU*ESC IM» 

All!  SHEAR? 

27fc 

PE*Ph«  ( 1  ,-ESC  (N»  9S*(DH/ESC(N.ll-l.)/2.IOHnsC(N,  9I*EH  SHEAR? 

?77 

_ 

PaAMAXllPEiO.) 

SHEAR? 

?7R 

sx«o. 

SHFAR? 

279 

ST«0. 

SHEAR? 

?R0 

...... 

S2«0. 

SHFAR? 

?fl  1 

SXTaO. 

SHEAR? 

282 

I H 3a 2 5 

SHEAR? 

?83 

IF  (NCALL  .6E.  91  00. TO  900 

SHEAR? 

?R* 

RETURN 

SHFAR? 

?85 

c  •  • 

•  •  • 

SHEAR? 

?HR 

c 

SHEAR? 

?8  l 

c 

PRINT  OARAGE  ARRAYS 

SHEAR? 

?RE 

c 

SHEAR? 

?B9 

C  •  • 

•  •  • 

SHE  AR? 

?90 

900 

IF  (EN  .CO.  o.)  RETURN 

SHtAR? 

?91 

PRINT  8DOO*K,J,lH3.ROTtENtTAUtEP 

SHFAR? 

?9? 

JNaO 

SHEAR? 

293 

00  1000  NQal'NANQ 

SHFAR? 

?9* 

NSaNSIZt (M«N0! #2 

SHFAR? 

?95 

IF  (NS  .EO.  0  .OR.  CN| JN*1 1  .£0.  O.ISO  TO  1100  SHEA*? 

?96 

IIaJN.NS-1 

SHFAR? 

297 

CNA(II)aCN(Hl 

SHEAR? 

296 

_ 

IF  (NS  .LE.  ?)  00  T0.9JS 

SHFAR? 

299 

00  970  la9(NS.2 

SHEAR? 

300 

I IaNS-1* JN 

SHEAR? 

301 

970 

CNA (It) aCN ( I I ) «CNA (11*21 

SHEAR? 

30? 

975 

CONTINUE 

SHE  AH? 

3d  3 

- 

PRINT  0500 tNO 

SHFAR? 

30* 

PRINT  9U01«(CNA(UN»Il»Ial(NSf2l 

SHEAR? 

305 

PRINT  9002* (CN| JN*II *Ia?*NS*2> 

SHEAR? 

30R 

•000 

FORMAT  OH  KalStSH  JaJStSM  |M9alS*5M  R0T*E16 

•3*AH  ENbEIO.9*  SHEAR? 

307 

: 

1  SH  TAuaE10.3*4H  EP-EU.3) 

SHEAR? 

308 
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FIGURE  42  LISTING  OF  SUBROUTINE  8HEAR2  (Conducted) 
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